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We revisit the problem of interacting electrons hopping on a two-leg ladder. A perturbative 
renormalization group analysis reveals that at half-filling the model scales onto an exactly soluble 
Gross-Neveu model for arbitrary finite-ranged interactions, provided they are sufficiently weak. The 
Gross-Neveu model has an enormous global 5*0(8) symmetry, manifest in terms of eight real Fermion 
fields which, however, are highly non-local in terms of the electron operators. For generic repulsive 
interactions, the two-leg ladder exhibits a Mott insulating phase at half-filling with d-wave pairing 
correlations. Integrability of the Gross-Neveu model is employed to extract the exact energies, 
degeneracies and quantum numbers of all the low energy excited states, which fall into degenerate 
SO{8) multiplets. One 5*0(8) vector includes two charged Cooper pair excitations, a neutral s = 1 
triplet of magnons, and three other neutral s = particle-hole excitations. A triality symmetry 
relates these eight two-particle excitations to two other degenerate octets which are comprised of 
sing/e-electron like excitations. In addition to these 24 degenerate "particle" states costing an energy 
(mass) m to create, there is a 28 dimensional antisymmetric tensor multiplet of "bound" states with 
energy y/Sm. Doping away from half-filling liberates the Cooper pairs leading to quasi-long-range 
d-wave pair field correlations, but maintaining a gap to spin and single-electron excitations. For very 
low doping levels, integrability allows one to extract exact values for these energy gaps. Enlarging 
the space of interactions to include attractive interactions reveals that there are four robust phases 
possible for the weak coupling two-leg ladder. While each of the four phases has a (different) 50(8) 
symmetry, they are shown to all share a common 50(5) symmetry - the one recently proposed by 
Zhang as a unifying feature of magnetism and superconductivity in the cuprates. 



I. INTRODUCTION 

Since the discovery of the cuprate superconductorsEI 
there has been renewed interest in the behavior of weakly 
doped Mott insulators There are two broad classes 
of Mott insulators, distinguished by the presence or ab- 
sence of magnetic order. More commonly spin rota- 
tional invariance is spontaneously broken, and long-ranaa 
magnetic order, typically antiferromagnetic, is realized.El 
There are then low energy spin excitations, the spin 1 
magnons. Alternatively, in a spin- liquid Mott insulator 
there are no broken symmetries, the magnetic order is 
short-ranged and there is a gap to all spin excitations : 
a spin-gap. 

In the cuprates the Mott insulator is antiferromagnet- 
ically ordered, but upon doping with holes the antifer- 
romagnetism is rapidly destroyed, and above a certain 
level superconductivity occurs. Below optimal doping 
levels, there are experimental signs of a spin gap opening 
at temperatures welL above the transition into the super- 
conducting phascOiS The apparent connection between a 
spin-gap and superconductivity has been a source of mo- 
tivation to search for Mott insulators of the spin-liquid 
variety. 

Although spin-liquids are notoriously difficult to 
achieve in two-dimensions,l3 it was realized that quasi- 
one-dimensional ladders would be more promising. .Par- 
ticular attention has focussed on the two-leg ladder .E3 At 
half-filling in the Mott insulator, the spin excitations can 
be described by a Heisenberg antiferromagnet, and due to 



the tendency for singlet bond formation across-jlJi&-jj4«igs 
of the ladder, spin-liquid behavior is expected Jj'oOlla In 
the past several years there have been extensiye,ranaly- 
ses of two-leaJajdders, particularly the Hubbardll^llil and 
t-J models,0il3ilZI both at half-filling and with doping. 
Based on numerical methods, including. Monte Carlo and 
density matrix renormalization groupt^a as well as ana- 
lytic approaches at weak coupling,E3Ej the basic behav- 
ior is established. At half-filling there is a spin-liquid 
phase with a spin-gap. Upon doping, the spin-gap sur- 
vives, although smaller in magnitude, and the system 
exhibits quasi-long-range superconducting pairing corre- 
lations, with approximate d-wave symmetry. This behav- 
ior is reminiscent of that seen in the underdoped cuprate 
superconductors. 

There are a number of experimental systems which can 
be described in terms of coupled two- leg ladders, "VjiliA 
exhibit a spin-gap in the insulating compound.LJL^I 
These materials are often very difficult to dope. In one 
case, doping has apparently been achieved, and under 
a pressuxe of 3GPa supercoruiuctivity is observed below 
12K.E£|lEI Carbon nanotubesEil constitute another novel 
material, wiiich can be modelled in terms of a two-leg 
ladder.Ej LJ Specifically, the low energy electronic exci- 
tations propagating down a single-walled nanotube can 
be mapped onto a two-leg ladder model with very weak 
interactions, inversely proportional to the tube radius. 

An obvious advantage of such low-dimensional corre- 
lated electron systems is (relative) theoretical simplic- 
ity. Indeed, in one-dimension many correlated elec- 
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tron models, including the Hubbard model, are exactly 
soluble.E3 Unfortunately, the Mott insulating phases of 
these one-dimensional models typically have gapless spin- 
excitations, and upon doping do not exhibit pairing. To 
date, we are unaware of any exactly soluble two-leg lad- 
der models which exhibit a gapped spin-liquid ground 
state. 

In this paper, we revisit models of interacting elec- 
tron hopping on a two-leg ladder, focusing on the behav- 
ior near half-filling. For generic short-range potentials, 
we derive a perturbative renormalization group valid for, 
weak interactions, much smaller than the bandwidth.E3't£l 
Remarkably, at half-filling the renormalization group 
transformation scales the system towards a special model 
with eMH-mous symmetry - the SO{8) Gross-Neveu (GN) 
model.E^ Scaling onto the GN model occurs independent 
of the initial interaction parameters, provided they are 
weak and predominantly repulsive. Thus, for weakly in- 
teracting two-leg ladders at half-filling universal low en- 
ergy properties are expected. Specifically, all properties 
on energy scales of order a characteristic GN mass (gap) 
m and distance scales longer than or of order v /m (where 
V is the Fermi velocity) are universal and determined by 
the GN model. In terms of microscopic parameters, the 
GN mass is of order m ~ te~'/^, where t is the Id band- 
width and U is a typical interaction strength, but is more 
profitably treated, along with v, as a phenomenological 
parameter. The universality predicted by the renormal- 
ization group can be profitably p^xploited because the 
S0{8) GN model is integrablejEr'B so that many of these 
universal properties can be computed exactly. To our 
knowledge, this is the first integrable model for a Mott- 
insulating spin liquid. It describes a state we call the D- 
Mott phase, because the Mott-insulator has short-range 
pairing correlations with approximate d-wave symmetry. 
We now summarize the results obtained from the 50(8) 
GN field theory. 

The primary inputj-frpm integrability is the complete 
excitation spectrum.EZrE3 The excitations of the GN 
model arc comprised of "particles" (i.e. sharp excita- 
tions with a single-valued energy-momentum relation) 
organized into 5*0(8) multiplets, as well as continuum 
scattering states of these particles. As expected for a 
Mott-insulating spin- liquid with no broken symmetries, 
each of these excitations is separated from the ground 
state by a non-zero gap. The lowest-lying particles 
come in three octets, all with mass m, i.e. dispersing 
as ei(g) = ^JrrP^^vq^ , where q is the deviation of the 
particle's momentum from its minimum energy value. 
One vector multiplet (conveniently denoted formally by 
a vector of Major ana fermions 77^, A = 1...8) con- 
sists entirely of collective two-particle excitations: two 
charge ±2e "Cooper pairs" around zero momentum, a 
triplet of spin-one "magnons" around momentum (tTjTt), 
and three neutral spin-zero "charge-density-wave" (or 
particle-hole pair) excitations. 50(8) transformations 
rotate the components of the vector into one another, 



unifying the pair, magnon, and charge-density-wave ex- 
citations. Indeed, the 50(5) subgroup rotating only the 
first five components of this vectopJs exactly the sym- 
metry proposed recently by Zhangiil to unify antiferro- 
magnetism and superconductivity in the cuprates. This 
vector octet, referred to as "fundamental" fermions in the 
field-theopt, literature, is related by a remarkable triality 
symmetryc3^c3 (present in the SO{N) GN model only for 
iV 8) to two other mass m octets: spinor and isospinor 
multiplets, called the even and odd kinks. These sixteen 
particles have the quantum numbers of individual quasi- 
electrons and quasi-holes. The triality symmetry thus 
goes beyond the 50(8) algebra to relate single-paijtij(jla 
and two-particle properties in a fundamental waycSH 
This relation also implies that pairing is present even in 
the Mott-insulator: the minimum energy to add a pair of 
electrons (as a member of the 50(8) vector multiplet) is 
m, reduced by a binding energy of m from the cost of 2m 
needed to add two quasi-electrons far apart. At energies 
above the 24 mass m states, there exists an antisymmet- 
ric tensor multiplet of 28 particles with mass -s/Sm. Each 
can be viewed as bound states of two different funda- 
mental fermions (or equivalently, two even or two odd 
kinks) . In this way their quantum numbers can be easily 
deduced by simple addition. The tensor states contribute 
additional sharp (delta-function) peaks to various spec- 
tral functions, providing, for instance, the continuation 
of the magnon branch near momentum (0, 0). For conve- 
nience, the quantum numbers (charge, spin, and momen- 
tum) of the vector and tensor excitations are tabulated 
in Tables 1 and 2. Finally, continuum scattering states 
enter the spectrum above the energy 2m. 

Combining the excitation spectrum of the GN model 
with the non-interacting spectrum and some additional 
arguments, we have also constructed schematic forms for 
several correlation functions of interest. In particular, 
in Sec. ^ we give detailed predictions and plots of the 
single particle spectral function (measurable by photoe- 
mission), the spin spectral function (measurable by in- 
elastic neutron scattering), and the optical conductivity. 
Integrability implies, for instance, sharp magnon peaks 
in the spin structure factor at k = (tTjTt), (0,0), and 
(±(fcj7i — kp2),T^) with minimum energy m, \/3m and 
^/Sm respectively (here kpi and kp2 are the Fermi mo- 
menta of the non-interacting system). Complete details 
can be found in Sec. ^ The optical conductivity has 
three principal features: a Drude peak around zero fre- 
quency, with exponentially small weight (~ e^™/"^) at 
low temperature; an "exciton" peak around uj = \/3m, 
exponentially narrow at low temperatures; and a contin- 
uum for Lu > 2to, due to unbound quasi-particle quasi- 
hole pairs. See Sec. g for more details and a figure. 

Our next calculations concern the relation of these re- 
sults to a recent study of microscopically 50(5) invariant 
ladder models by Scalapino, Zhang and Hanke (SZH).lj 
These authors consider the strong coupling limit of a 
certain locally-interacting two-leg ladder model designed 
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to exhibit exact 5*0(5) symmetry. Their model has an 
on-site interaction \U\ S> an intra-rung interaction 
\V\ 3> i, and a magnetic rung-exchange interaction J, 
related to one another by the 5*0(5) symmetry. In the 
U-V plane they derive a strong-coupling phase diagram, 
including the case of attractive interactions with U and 
V negative. We have analyzed general 50(5) invariant 
two-leg ladder models in the opposite limit of weak in- 
teractions, deriving as a special case the corresponding 
weak-coupling phase diagram for their model. In fact, al- 
though we have not explored the full 9-dimcnsional space 
completely, for all bare couplings we have considered, in- 
cluding attractive interactions that break 50(5) symme- 
try explicitly, the RG scales the system into the 50(5) 
subspace. When the interactions are predominantly re- 
pulsive, the 5*0(5) system falls into the basin of attrac- 
tion of the D-Mott phase, and the above results apply. 
As negative interactions are introduced, four other phases 
emerge: an S-Mott spin-liquid, with short-range approx- 
imate s-wave pairing symmetry, a charge-density-wave 
(CDW) state with long-range positional order at (tTjTt), 
a spin-Peierls phase with kinetic energy modulated at 
(tt, tt), and a Luttinger liquid (C2S2, in the nomencla- 
ture of Ref. |l^) phase continuously connected to the non- 
interacting system. The first two of these also occur in 
the strong-coupling limit, though their positions in the 
phase diagram (Fig. 10) are modified. The phase dia- 
grams at weak and strong coupling differ in non-trivial 
ways, implying a rather complex evolution of the system 
with increasing U and V. In weak coupling, all four non- 
trivial phases have distinct asymptotic 50(8) symme- 
tries, enhanced from the common bare 5*0(5). Further- 
more, critical points describing the transitions between 
the various phases can also be identified. In particular, 
the D-Mott to S-Mott and CDW to spin-Peierls criti- 
cal points are c = 1 conformal field theories (single mode 
Luttinger liquids) , which in weak-coupling are acccompa- 
nied by a decoupled massive 50(6) sector. The S-Mott 
to CDW and D-Mott to spin-Peierls transitions are Ising 
critical theories (c = 1/2), with decoupled massive 50(7) 
sectors in weak-coupling. There is also a multi-critical 
point describing a direct transition from the D-Mott to 
CDW or from the S-Mott to spin-Peierls phases, which is 
simply a product of the c — 1 and c = 1/2 critical points. 

Our final results concern the effects of doping a small 
density of holes (or electrons) into the D-Mott spin- liquid 
phase at half-filling. For very small hole concentrations, 
the modifications of the Fermi velocities by band curva- 
ture effects can be ignored, and the doping incorporated 
simply by including a chemical potential term coupled 
to the total charge Q in the GN model; iJ^ — K.— jiQ. 
An analogous procedure is employed by Zhan^iJ in his 
study of the 5*0(5) non-linear sigma model. Because the 
charge Q is a global 50(8) generator, integrability of the 
GN model is preserved, and furthermore many of the 
50(8) quantum numbers can still be employed to label 
the states. We find that doping occurs only for 2/i > to, 
at which point Cooper pair "fundamental fermions" en- 



ter the system and effectively form a Luttinger liquid 
with a single gapless charge mode (with central charge 
c = 1). This phase (often denoted "ClSOjli jsiiU has a 
gap to spin excitations. Previous worklli'EjO has ap- 
proached this phase via controlled perturbative calcu- 
lations in the interaction strength, at fixed doping x 
away from half-filling. Here, we are considering a dif- 
ferent order of limits, with fixed (albeit weak) interac- 
tions in the small doping limit, a; — s- 0. In this limit, the 
Cooper-pairs being dilute behave as hard-core bosons or 
free fermions. Although the spin-gap is preserved in the 
doped state, it is discontinuous as a; — > 0"*". The discon- 
tinuity can be understood as the binding of an inserted 
spin-one magnon to a Cooper pair in the system to form 
a mass \/3m tensor particle, reduced by the binding en- 
ergy (2 — \/3)m from its bare energy. The spin-gap thus 
jumps from As{x ^ 0) ~ m to As{x = 0+) = — l)m 
upon doping. Such binding of a pair to a magnon has 
been observed numerically in tath Hubbard and t-J lad- 
ders by Scalapino and White .e3 Similarly, the energy to 
add an electron (for the hole-doped system) jumps from 
Ai_(a: = 0) = 3to/2 to Ai_(a; = 0+) = to/2, the same 
as the energy to add a single hole. When many pairs are 
present, we have not succeeded in obtaining exact ex- 
pressions for the spin and single-particle gaps, but argue 
that the spin gap should decrease with increasing dop- 
ing, since the added magnon is attracted to an increasing 
density of Cooper pairs. It seems likely, however, that in- 
tegrability could be exploited even in this case to obtain 
exact results, and hope that some experts may explore 
this possibility in the future. 

Finally, we briefly address the behavior of the spin- 
spectal function for the doped ladder atj^ergies above 
the spin-gap. In a recent preprint SZHa have argued 
that in this regime the spin-spectral function for a model 
with exact SO (5) symmetry should exhibit a sharp res- 
onance at energy 2/i and momentum (tt, tt), the so- 
called TT— resonance (introduced originally by Zhang to 
explain the 42meV neutron scattering peak in the super- 
conducting Cuprates). We show that a delta- function 
TT— resonance requires, in addition to SO (5) symmetry, 
the existence of a non-zero condensate density in the 
superconducting phase. Since condensation is not pos- 
sible in one-dimension, this precludes a delta-function. 
TT— resonance. Following a recent suggestion by Zhang,c3 
we address briefly the possibility of a weaker algebraic 
singularity in the spin spectral function. Regardless of 
the nature of the behavior in the vicinity of w = 2fi, we 
expect spectral weight at energies below 2/i but above 
the spin-gap discussed above. 

The remainder of the paper is organized as follows. In 
Sec. H we describe the model Hamiltonian for the inter- 
acting ladder, reduce it to the continuum limit, bosonize 
the nine distinct interaction channels, and apply the 
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renormalization group (RG) transformation. Sec 
details the simplifications that occur upon RG scaling, 
presents the bosonized form of the Hamiltonian in the 
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D-Mott phase, and for completeness demonstrates the 
short-range d-wave correlations found in Ref . |l^ The bulk 
of the field-theoretic analysis is contained in Sec. [V, By 
refermionizing the bosonized hamiltonian, we obtain the 
GN model exposing the exact SO (8) symmetry, and de- 
scribe why this symmetry is hidden in the original vari- 
ables. The triality symmetry is identified, and used to 
understand the degeneracy between the three mass m 
octets. To help in developing an intuition for the GN 
model, several approximate pictures are presented to un- 
derstand the excitations: a mean field theory which is 
asymptotically exact for iV — » oo in a generalized SO{N) 
GN model, and a semi-classical theory based on the 
bosonized (sine-Gordon-like) form of the Hamiltonian. 
We conclude Sec. |^ by proving the uniqueness of the 
ground-state in the D-Mott phase and determining the 
quantum numbers of the 24 -1- 28 = 52 particles. The 
latter task is complicated by the necessity of introducing 
Jordan- Wigner strings, which are required to preserve 
gauge-invariance under an unphysical gauge symmetry 
introduced in bosonization. The string operators mod- 
ify the momenta of the certain excitations by a shift of 
(tt, tt) from their naive values determined from the GN 
fermion operators. With the field-theoretic analysis com- 
plete, we go on to discuss correlation functions in Sec. 0, 
giving detailed predictions for the single-particle spectral 
function, spin spectral function, optical conductivity, and 
various equal-time spatial correlators. Sec. Vl describes 



the construction of general 5*0(5) invariant models in 
weak-coupling, their phases, and the phase-diagram of 
the Scala pino- Zhang-Hanke model in weak-coupling. Fi- 
nally, Sec. VII describes the behavior of the D-Mott phase 
upon doping, including the behavior of various gaps, and 
a discussion of the status of the SO (5) "tt resonance" in 
one dimension. Various technical points and long equa- 
tions are placed in the appendices. Appendix A gives the 
full set of nine RG equations at half-filling. Appendix B 
discusses gauge redundancy and the multiplicity of the 
ground state in different phases. Appendix C constructs 
spinor and vector representations of 5*0(5), Appendix D 
relates 50(5) and 50(8) currents, and Appendix E gives 
the five RG equations in the reduced 5*0(5) subspace. 



II. MODEL 

We consider electrons hopping on a two-leg ladder as 
shown in Fig. 1. In the absence of interactions, the 
Hamiltonian consists of the kinetic energy, which we as- 
sume contains only near-neighbor hopping. 
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ta\^{x + l)a^Jx) + {1 ^ 2) 



-tj_ala{x)a2a{x) + h.c. 



(2.1) 



where af(a|) is an electron annihilation (creation) oper- 
ator on leg £ of the ladder (i' = 1, 2), a; is a discrete co- 



ordinate running along the ladder and a =t, i is a spin 
index. The parameters t and t± are hopping amplitudes 
along and between the leg's of the ladder. 

Being interested in weak interactions, we first diago- 
nalize the kinetic energy in terms of bonding and anti- 
bonding operators: = (0,1,0 + {— 0,2,0) /V^, with 
i = 1,2. The Hamiltonian is then diagonalized in mo- 
mentum space along the ladder, describing two decou- 
pled (bonding and anti-bonding) bands. Focussing on 
the case at half-filling with one electron per site, both 
bands intersect the Fermi energy (at zero energy) pro- 
vided t± < 2t. Moreover, due to a particle/hole symme- 
try present with near neighbor hopping only, the Fermi 
velocity Vi in each band is the same, denoted hereafter as 
V. It is convenient to linearize the spectrum around the 
Fermi points at ±kpi (see Fig. 1), which at half-filling 
satisfy kpi + fcjr2 = Upon expanding the electron op- 
erators as. 
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(2.2) 



the effective low energy expression for the kinetic energy 
takes the form, Hq — J dxTLo, with Hamiltonian density. 



(2.3) 



This Hamiltonian describes Dirac Fermions, with four 
flavors labelled by band and spin indices. Since all 
flavors propagate both to the right and left with the 
same velocity, the model exhibits an enlarged symme- 
try. Specifically, if the four right (and left) moving Dirac 
Fermions are decomposed into real and imaginary parts, 
V-P^a = i^Uo + ^C^o)/V2 where P = R/L and are 
Majorana fields, the eight right (and left) moving Ma- 
jorana fields, denoted ^p^i with A = 1,2,..., 8 form an 
eight component vector. The Hamiltonian density, when 
re-expressed in terms of these eight component vectors 
takes the simple form 



V 



8 



(2.4) 



which is invariant under independent global 50(8) rota- 
tions among either the right or left vector of Majorana 
fields. This enlarged 0(8)^ x 0(8)l symmetry is only 
present at half- filling with particle/hole symmetry. 
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FIG. 1: A two-leg ladder and its band structure. In the 
low-energy limit, the energy dispersion is linearized near 
the Fermi points. The two resulting relativistic Dirac 
Fcrmions are distinguished by pseudospin indices i — 1,2 
for the anti-bonding and bonding bands, respectively. 

Electron-electron interactions scatter right-moving 
electrons into left-moving electrons and vice-versa, de- 
stroying this large symmetry. For general spin- 
independent interactions the symmetry will be broken 
down to U{1) X SU (2), corresponding to total charge and 
spin conservation. In the following we consider general 
finite-ranged spin-independent interactions between the 
electrons hopping on the two-leg ladder. We assume the 
typical interaction strength, [/, is weak - much smaller 
than the bandwidth. We focus on the effects of the inter- 
actions to leading non-vanishing order in U. In this limit 
it is legitimate to keep only those pieces of the interac- 
tions which scatter the low energy Dirac Fermions. Of 
these, only those involving four-Fermions are marginal, 
the rest scaling rapidly to zero under renormalization. 
Moreover, four-Fermion interactions which are chiral, say 
only scattering right movers, only renormalize Fermi ve- 
lo(jit|iea and can be neglected at leading order in small 
f/.EallB All of the remaining four-Fermion interactions can 
be conveniently expressed in terms of currents, defined as 



(2.5) 



lij = CiatapCjp, lij — - Cia{e(T)apCjp, (2.6) 

where the i?, L subscript has been suppressed. Both J 
and / are invariant under global SU{2) spin rotations, 
whereas J and I rotate as SU{2) vectors. Due to Fermi 
statistics, some of the currents are (anti-)symmetrical 



(2.7) 



so that In — (no sum on i). 

The full set of marginal momentum-conserving four- 
Fermion interactions can be written 



-w(i) - hP 7 ■ T ■ 



h'^ J ... J . . 

^fij^RiiJ^jj — fijJnii ' JLjj- 



(2.8) 



Here fij and bij denote the forward and backward 
(Cooper) scattering amplitudes, respectively, between 
bands i and j. Summation on i,j = 1,2 is implied. 
To avoid double counting, we set fa = (no sum on 
i). Hermiticity implies 612 = &21 and parity symmetry 
{R ^ L) gives /12 = /21, so that there are generally 
eight independent couplings h^'-^ , 622^, 6^2'^, and fi2- At 
half- filling with particle/hole symmetry 611 = &22- Addi- 
tional momentum non-conserving Umklapp interactions 
of the form 
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u''-!^ -I - 

^IJ RIJ LIJ 



U^ I^ - - 

^IJ RtJ 



LIJ 



h.c. 



(2.9) 



are also allowed, (here 1 = 2, 2 = 1). Because the 
currents (lij), lij are (anti-) symmetric, one can always 
choose U12 = U21 for convenience. We also take = 
since lu = 0. With particle/hole symmetry there are 
thus just three independent Umklapp vertices, u^i, u^j, 
and U12. Together with the six forward and backward 
vertices, nine independent couplings are required to de- 
scribe the most general set of marginal non-chiral four- 
Fermion interactions for a two-leg ladder with parti- 
cle/hole symmetry at half-filling. 

Since our analy sis below makes heavy use of abelian 
Bosonization,E3E3 it is convenient at this stage to con- 
sider the Bosonized form of the general interacting the- 
ory. To this end, the Dirac fermion fields are expressed 
in terms of Boson fields as 



(2.10) 



where P = R/ L = ±. To ensure that the Fermionic oper- 
ators anti-commute the Boson fields are taken to satisfy 

[(l)pia{x), (I)pjf3{x')] = iPTT5i-j5apSgn{x - x'), (2.11) 
[(l>Riaix),(j)Ljl3ix')] = inSij6ai3 . (2.12) 



Klein factors, satisfying 



(2.13) 



have been introduced so that the Fermionic operators in 
different bands or with different spins anticommute with 
one another. 

It will also be convenient to define a pair of conjugate 
non-chiral Boson fields for each flavor. 



V^ZQ — 4'Ria ^" ^Lia 7 
^ia = ^Ria ^Lia-t 

which satisfy 

[tpix),e{x')] = -iinOix' - x). 



(2.14) 
(2.15) 



(2.16) 



Here, and in the remainder of the paper, we denote by 
Q{x) the heavyside step function to avoid confusion with 
the 9 fields defined in Eq. 2.15 above. The field 9ia is a 



displacement (or phonon) field and ipia is a phase field. 

is 



The Bosonized form for the kinetic energy Eq. 2.3 



on ^ — ' 



(2.17) 



which describes density waves propagating in band i and 
with spin a. 

This expression can be conveniently separated into 
charge and spin modes, by defining 



(2.18) 
(2.19) 
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and similarly for ip. The -\/2 ensures that these new fields 
satisfy the same commutators, Eq. ( 2.16 ). It is also con- 
venient to combine the fields in the two bands into a ± 
combination, by defining 



{din ± i 



2, 



(2.20) 



where fi — p,a, and similarly for tp. It will sometimes be 
convenient to employ charge/spin and flavor decoupled 
chiral fields, defined as 



(2.21) 



with P = R/L = ±. 

The Hamiltonian density Tio can now be re-expressed 
in a charge/spin and flavor decoupled form, 

Ho = ^ y^me^if + (a,^^±)2]. (2.22) 

The fields (?p+ and ipp^ describe the total charge and cur- 
rent fluctuations, since under Bosonization, cl,^^Cp^^ = 
d^9p+/Tr and vPcl^^c^i^ = d^^pp+lir. 

The interaction Hamiltonians can also be readily ex- 
pressed in terms of the Boson fields. The momentum 
conserving terms in Eq. 2.8 can be decomposed into two 
contributions, = H}^ + Ti.\^''\ the first 2 involving 
gradients of the Boson fields. 



^(la) ^ 



1 



167r2 ^ 



^ A^^iid^e^^f - id,^^±)% (2.23) 



with coefficient Ap± — 2(cjj ± f^2) ^^^^ Aa± — — (c^^ ± 
/f2)/2, whereas the second contribution involves cosines 
of the Boson fields: 

n'-j-''^ = -2r5?2 cos ipp^ cos 6*^+ 

+ cos e„+ (26^1 cos 9^^ + 2r/f2 cos (p^- ) 

— cos (/3p_ (r6]'2 cos 0cr- + &12 ^OS (^cr- ) , (2.24) 

with — ^12 i ^^^12 ■ Similarly, the Umklapp interac- 
tions can be Bosonized as. 



— 16Tu1^ cos 0p-)_ ccj 
cos 6p+ {2u^2 '^'^^ d<y 



cos ipp- 



4u^2 cos Op+ cos 6a+ 
2ruj;2COS(^^_), (2.25) 

with = ± 4u^2- Here T — Ki|Ki|K2|K2i is a prod- 
uct of Klein factors. Since = 1, we can take F = ±1. 
Hereafter, we will put F = 1. 

In the absence of electron-electron interactions, the 
Hamiltonian is invariant under spatially constant shifts 
of any of the eight non-chiral Boson fields, 6p± and ^Pp±. 
With interactions jive of the eight Boson fields enter 
as arguments of cosines, but for the remaining three - 
(pp+ , and Op- - this continuous shift symmetry is 
still present. For the first two fields, the conservation law 
responsible for this symmetry is readily apparent. Specif- 
ically, the operators exp(iaQ) and exp(mS'z), with Q the 



total electric charge and Sz the total z-component of spin, 
generate "translations" proportional to a in the two fields 
(Pp+ and ^Pa+- To see this, we note that Q — J dxp{x) 
with p{x) — dxOp+li^ the momentum conjugate to V'p+i 
whereas Sz can be expressed as an integral of the mo- 
mentum conjugate to ^Pa+- Since the total charge is con- 
served, [Q,H] = 0, the full Hamiltonian must therefore 
be invariant under ipp^ — > Lpp+ + a for arbitrary constant 
a, precluding a cosine term for this field. Similarly, con- 
servation of Sz implies invariance under ^pa+ 'Pa+ + o- 
The conservation law responsible for the symmetry un- 
der shifts of the third field, 9a-, is present only in the 
weak coupling limit. To see this, consider the opera- 
tor, V = kpiJi + kp2J2, with J., = Y^ai^aia " N ^la) , 

where Npia is the total number of electrons in band i 
with spin a and chirality P. At weak coupling with 
Fermi fields restricted to the vicinity of fc^^, this oper- 
ator is essentially the total momentum. Since the total 
momentum is conserved up to multiples of 27r, one has 
AV — ±27rrt = ±2n(A;pi -f- kp2) for integer n. Moreover, 
since the Fermi momenta kpi are in general unequal and 
incommensurate, this implies that A Ji = A J2 = ±2n, or 
equivalently that Ji — J2 is conserved at weak coupling. 
Since Ji — J2 = / dxj(x) with j(x) = d^'^p-j'K the mo- 
mentum conjugate to dp-, this conservation law implies 
invariance under dp- 6p- + a. 

The remaining 5 Boson fields, entering as arguments 
of various cosine terms, will tend to be pinned at the 
minima of these potentials. Two of these 5 fields, 9a- 
and (fa- , are dual to one another so that the uncertainty 
principle precludes pinning both fields. Since there are 
various competing terms in the potential seen by these 
5 fields, minimization for a given set of bare interaction 
strengths is generally complicated. For this reason we 
employ the weak coupling perturbative renoipaalization 
group transformation, derived in earlier workJl3'E3 Upon 
systematically integrating out high-energy modes away 
from the Fermi points and then rescaling the spatial co- 
ordinate and Fermi fields, a set of renormalization group 
(RG) transformations can be derived for the interaction 
strengths. Denoting the nine interaction strengths as gi, 
the leading order RG flow equations take the general 
form, digi = Aijkgjgk, valid up to order g^. For com- 
pleteness the RG flow equations are given explicitly in 
Appendix A. Our approach is to integrate the RG flow 
equations, numerically if necessary, to determine which 
of the nine coupling constants are growing large. 

Under a numerical integration of these nine flow equa- 
tions it is found that some of the couplings remain small, 
while others tend to increase, sometimes after a sign 
change, and then eventually diverge. Quite surprisingly, 
though, the ratios of the growing couplings tend to ap- 
proach fixed constants, which are indepependent of the 
initial coupling strengths, at least over a wide range in 
the nine dimensional parameter space. These constants 
can be determined by inserting the Ansatz, 
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9i 



(2.26) 



into the RG flow equations, to obtain nine algebraic equa- 
tions quadratic in the constants gio- There are various 
distinct solutions of these algebraic equations, or rays in 
the nine-dimensional space, which correspond to differ- 
ent possible phases. But for generic repulsive interactions 
between the electrons on the two-leg ladder, a numerical 
integration reveals that the flows are essentially always 
attracted to one particular ray. In the next sections we 
shall consider the properties of this phase, which, for rea- 
sons which will become apparent, we denote by D-Mott. 



III. D-MOTT PHASE 

In the phase of interest, two of the nine coupling con- 
stants, b'^^ and /fj, remain small, while the other seven 
grow large with fixed ratios: 



K2 = J&?2 



fP 
J 12 



= .9 > 0. 



(3.1) 
(3.2) 



Once the ratio's are fixed, there is a single remaining cou- 
pling contant, denoted g, which measures the distance 
from the origin along a very special direction (or "ray" ) 
in the nine dimensional space of couplings. The RG equa- 
tions reveal that as the flows scale towards strong cou- 
pling, they are attracted to this special direction. If the 
initial bare interaction parameters are sufficiently weak, 
the RG flows have sufficient "time" to renormalize onto 
this special "ray" , before scaling out of the regime of per- 
turbative validity. In this case, the low energy physics, on 
the scale of energy gaps which open in the spectrum, is 
universal, depending only on the properties of the physics 
along this special ray, and independent of the precise val- 
ues of the bare interaction strengths. 

To expose this universal weak coupling physics, we use 
Eq. to replace the nine independent coupling con- 
stants in the most general Hamiltonian with the single 
parameter g, measuring the distance along the special 
ray. Doing so reveals a remarkable symmetry, which is 
most readily exposed in terms of a new set of Boson fields, 
defined by. 



{9,^)2 = (0, ¥>).+ , 



(3.3) 



The first three are simply the charge/spin and flavor 
fields defined earlier. However, in the fourth pair of fields, 
6 and cp have been interchanged. It will also be useful to 
consider chiral boson fields for this new set, defined in 
the usual way. 



with a = 1,..,4, and P = R/L = ± as before. The 
first t hree of thes e chi ral fields satisfy the commutators 
Eq. ( |2.11| ) and { ^12^ . But for the fourth field, since 
(/)p4 — P4)pp-, the second commutator is modified to 
[0i?.4, 4>la\ = -«7r. 

In terms of these new fields, the full interacting Hamil- 
tonian density along the special ray takes an exceedingly 
simple form: Ti — Ti-o + Tii, with 



Wo = £^5]P.0a)' + (9.^a)'], 



(3.5) 



1^ cos 9a cos db- 



(3.6) 



We now briefly discuss some of the general physical prop- 
erties which follow from this Hamiltonian. In the next 
sections we will explore in detail the symmetries present 
in the model, and the resulting implications. 

Ground state properties of the above Hamiltonian can 
be inferred by employing semi-classical considerations. 
Since the fields cpa enter quadratically, they can be in- 
tegrated out, leaving an effective action in terms of the 
four fields 9a- Since the single coupling constant g is 
marginally relevant and flowing off to strong coupling, 
these fields will be pinned in the minima of the cosine 
potentials. Specifically, there are two sets of semiclassical 
ground states with all 9a = 2na7r or all 9a — {2na + l)7r, 
where Ua are integers. Excitations will be separated from 
the ground state by a finite energy gap, since the fields 
are harmonically confined, and instanton excitations con- 
necting different minima are also costly in energy. 

Since both 9a± fields are pinned, so are the spin-fields 
in each band, 9ia- {i — 1,2). Since dx9ia is proportional 
to the z-component of spin in band i, a pinning of these 
fields implies that the spin in each band vanishes, and 
excitations with non-zero spin are expected to cost finite 
energy: the spin gap. This can equivalently be inter- 
preted as singlet pairing of electron pairs in each band. 
It is instructive to consider the pair field operator in band 



-^(vio+s 



(3.7) 



With 9icr ^ 0, f ip can be interpreted as the phase of the 
pair field in band i. The relative phase of the pair field 
in the two bands follows by considering the product 



(3.8) 



0™ = i^a + P0a)l2, 



(3.4) 



with r = Ki|Ki|K2|K2i = 1- Since 9^ = (pp- the rel- 
ative phase is also pinned by the cosine potential, with 
a sign change in the relative pair field, A^AJ < 0, cor- 
responding to a D-wave symmetry. Being at half-filling, 
the overall charge mode, 9p+ is also pinned - there is 
a charge gap - and the two-point pair field correlation 
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function falls off exponentially with separation. We re- 
fer to this phase as a "D-Mott" phase, having D-wave 
pairing correlations coincident with a charge gap. Upon 
doping the D-Mott phase away from half-filling, gapless 
charge fluctuations are expected in the (p+) sector, and 
power-law D-wave pairing correlations develop. 

It is worth noting that the fully gapped D-Mott phase 
has a very simple interpretation in the strong cou- 
pling limit. Two electrons across each of the rungs of 
the two-legged ladder form singlets, of the usual form 
I T: i) ^ I i: T)j where the two states refer to electrons on 
leg 1 or 2, respectively. This two-electron state can be 
re-written in the bonding anti-bonding basis, and takes 
the form, | t^, — ) — | — , t^), where the two states now re- 
fer to bonding and anti-bonding orbitals. This resembles 
a local Cooper pair, with a relative sign change between 
bonding and anti-bonding pairs: an approximate D-wave 
symmetry. 



IV. 50(8) GROSS-NEVEU MODEL 

As shown above, the bosonized effective Hamiltonian 
on energy scales of order the gap is exceptionally simple 
in the D-Mott phase. In this section, we show that this 
simplicity is indicative of a higher symmetry, and explore 
its ramifications upon the spectrum. 



A. Gross-Neveu Model 



An obvious symmetry of the bosonic action, Eqs. 



3.6 , is permutation of the fields 9a PabOt, where Pat is 
a permutation matrix. In fact, this is only a small subset 
of the true invariances of the model. As is often the case, 
abelian bosonization masks the full symmetry group. It 
can be brought out, however, by a refermionization pro- 
cedure. We define "fundamental" (Dirac) fermion oper- 
ators ippa with a — 1,2,3,4 via 



a = 1...3 



(4.1) 



and P — R, L — ±1, as before. The Klein factors are 
given by 



(4.2) 
(4.3) 



In the re-Fermionization of the fourth field we have cho- 
sen to include a minus sign for the left mover. This is 
convenient, due to the modified commutators between 
the left and right fields: [(^ ^4, ipL ^] — —in, in contrast to 
the "standard" form in Eq. 2.12| . 

In these variables, the effective Hamiltonian density 
becomes 



where tpa = {fpRajtl^La), and r is a vector of Pauli ma- 
trices acting in the R, L space. Here, summation over 
repeated indices, a — 1,2,. .,4 is implicit. It is remark- 
able that the Hamiltonian can be written locally in the 
"fundamental" fermion variables, which are themselves 
highly uon-locally related to the "bare" electron opera- 
tors. 

A further simplification arises upon changing to Majo- 
rana fields. 



(4.5) 



The Hamiltonian density then takes the manifestly in- 
variant form 

n = \vnAid..VHA - ^VlaiO^Vla + gGf'Gt^, (4.6) 



where the currents are 

and ^, S = 1 . . . 8. 

B. 50(8) Symmetry 



(4.7) 



Eq. 4.6 is the standard form for the S'0(8) Gross- 



Neveu niC|deL-,\phich has been intensively studied in the 
literature B3e3^e3'H We first discuss its manifest symme- 
try properties. 

The 28 currents Gp^ generate chiral 5*0(8) transfor- 
mations. For g = 0, Eq. 4^ has two independent sym- 
metries under separate rotations of the left- and right- 
moving fields. For g 7^ 0, however, only simultaneous 
rotations of both chiralities are allowed. More precisely, 
the unitary operators 



(4.8) 



(4.4) 



generate global orthogonal transformations of the Majo- 
rana fields, 

U\x)VpaU{x) = 0^s{x)Vpb, (4.9) 
where the orthogonal matrix 0{x) is given by 

Q^^^ ^ ^iXABTAB ^ (4.10) 

Here the Tab {A > B) are the 28 generators of 50(8) 
in the fundamental representation, wit h m atrix elements 
[Tab\co = ^i^Ac^BD - <5ad<5bc)/2. Eq. iudicatcs that 
the rjpA transform as S'0(8) vectors. Similarly, the cur- 
rents Gp^ are rank 2 S'0(8) tensors. 

It is worth noting that despite the non-local relation 
between the fundamental and bare fermion operators, the 
50(8) symmetry remains local in the bare electron basis. 
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This follows from the fact that the chiral SO{8) currents 
in the two bases are actually linearly related, i.e. 



(4.11) 



where G^^ = i£,PA£.PB, and the bare Majorana operators 
are defined by 



Cpil = (Cp4 + i^ps) , 
Cp2T = (^p6 + «Cp5) , 

Cp2l = (Cp8 + ■iS.pj) ■ 



(4.12) 
(4.13) 
(4.14) 
(4.15) 



The precise forms of the tensors Mp are complicated and 
not particularly enlightening. Nevertheless, the existence 
of the linear relation between currents implies that the 
unitary operator U{x) also generates local rotations of 
the bare electron fields. In these variables, however, the 
50(8) symmetry is hidden, because Mj^ ^ M^, which 
implies different rotations must be performed amongst 
right- and left-moving electron operators. 

Finally, it is instructive to see how the conserva- 
tion of total charge and spin, corresponding to a global 
[/(I) X SU (2) symmetry, is embedded in the larger 5*0(8) 
symmetry. To this end, consider the total electron charge 
operator, Q, which in terms of the low energy fields can 
be written, 



g = 2 



dx^ = 2 / dxiGf + Gf), (4.16) 



where V'pi is a fundamental Gross-Neveu fermion. The 
U{1) charge symmetry is thus seen to be equivalent to 
the 50(2) symmetry of rotations in the 1 — 2 plane of 
the eight-dimensional vector space. Similarly, the total 
spin operator 



dx [J nix) + Ji,(x)], 



(4.17) 



with J p{x) — Jpii{x), can be re-expressed in terms of 
50(8) generators by using. 



^abc^bc 



(4.18) 



with a, 6, c = 3, 4, 5 = x, y, z. Thus we see the equivalence 
between the SU{2) spin rotations and 50(3) rotations in 
the 3-dimensional sub-space 3 — 4 — 5 of the eight dimen- 
sional vector space. Rotations in the five-dimensional 
subspace 1 — 2 — 3 — 4 — 5, correspond to global 50(5) 
rotations which unify the charge and spin degrees of free- 
dom. 

In the absence of interactions in the Gross-Neveu 
model, all of the excitations including spin remain mass- 
less. In this case there is an independent SU{2) spin 



symmetry in the right and left moving sectors. The spin 
currents J p can then be shown to satisfy, 



[j;{x), 4{x')] = 6{x - x')te'"''r^{x) + 
p 

i—k6ab6'{x - a;'), 

ZTT 



(4.19) 
(4.20) 



with a,b,c = x, y, z and k = 2. This is referred to as an 
SU{2) current algebra at level (k) two. 

We conclude this subsection by answering a question 
which may have occurred to the alert reader: why is 
the symmetry of the model 50(8) rather than 0(8)? 
Based on Eq. 4.6, it would appear that any transfor- 
mation of the form tjpa — > OabVpb would leave the 
Hamiltonian invariant, including improper rotations with 
detO — —1. The presence of such improper rotations 
means 0(8) — 50(8) x Z2, since any orthogonal ma- 
trix can be factored into a product of matrix with de- 
terminant one and a particular (reflection) matrix, e.g. 
Cab = ^AB — 2Sjii5bi- We have already shown above 
that the 50(8) symmetry is physical - i.e. the symme- 
try generators act within the Hilbert space of the phys- 
ical electrons. It is straightforward to show that the Z2 
reflection is however, unphysical. To see this, imagine 
performing the Z2 reflection effected by C above, which 
takes rjpi — > —rjp^. Using the bosonization rules, this cor- 
responds to 9i — > —61 and (pi — > — Returning to the 
physical fields, one finds that the bare electron operators 
transform much more non-trivially: 



(4.21) 



As we shall show in Sec. IV. E. 3, a single GN fermion op- 
erator, such as i/'pi, is unphysical. The Z2 reflection thus 
takes a physical electron operator into an unphysical one, 
which implies that the symmetry cannot be effected by 
a unitary operator within the Hilbert space of the elec- 
trons. For this reason, the true symmetry group of the 
ladder model is 50(8). 

GN Fermions 

(IwO'particle) 



La 



I I, 



TRr. 



Rfl 



Lfl 



¥' 

I Rn 



Odd Kinks Even Kinks 

f single-partick) (single-particle) 

FIG. 2: Triality between GN fermions, even kinks and 
odd kinks. The 50(8) GN Hamiltonian is identical in 
terms of these three sets of fermionic operators. Opera- 
tors in the gray areas are physical and gauge independent 
(see Sec. IV. E), while the other fermion operators must 
be "dressed" by an appropriate Jordan- Wigner string to 
remain in the physical Hilbert space. 
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C. Triality 

Most of the above properties hold more generally for 
the SO{N) GN model, even for N ^ 8. However, the 
case iV = 8 is extremely special, and in fact possesses an 
additional triality symmetry not found for other N . A 
useful reference is Ref. |4^ , ^ . 

To expose the additional symmetry, we return to the 
sine-Gordon formulation. Essentially, the triality opera- 
tion trades the original basis {Oa} in the four-dimensional 
space of boson fields for cither one of two other orthog- 
onal bases. Explicitly, the two alternate choices are the 
even and odd fields 0a^° , where 



qe/o _ 

qe/o _ 
^3 — 

ft- 



O3 ± 04)/2, 
O3 ± 0i)/2. 



(4.22) 
(4.23) 
(4.24) 
(4.25) 



Here the upper and lower signs apply to the even and odd 
fields, respectively, and identical definitions hold for the 
dual (fit/" and chiral (^p(,° bosons. The bosonizcd Hamil- 



tonian in Eqs. 3.5, 3.6 is invariant under either change of 
variables, i.e. 



H[ea] = H[ei] = H[0:] 



(4.26) 



For each of these bases, an inequivalent refermionization 
is possible, analogous t o th e introduction of the funda- 
mental fermions in Eq. 4.1. In particular, the Hamilto- 



nian is unchanged in form when rewritten in terms of 
either the even or odd fcrmion operators, 



Vpa 



(4.27) 



It should be noted that the set of even and odd fermion 
operators contains all the bare electron fields. In partic- 
ular. 







= C1.1T, 


(4.28) 


ru2 = 


Cr2T 


V'l2 = Cl2T, 


(4.29) 




CH2i 




(4.30) 








(4.31) 


eight 


even 


and odd fields (^^^ 


and '0^£j) are 



not simply related, however, to the electron fields. 



D. Conventional Gross-Neveu Excitation spectrum 

The SO{N) GN model is integrable, and the excitation 
spectrum is known exactly. To organize the presentation, 
we divide the discussion of the excitation spectrum into 
two parts. In this subsection, we summarize known re- 
sults for the conventional GN model. The precise nature 
of the excitations for the two-leg ladder model, however, 



differs from those in the conventional GN model. This 
difference arises from the non-local relation between the 
electron and GN fields. Excitations within the GN model 
must be slightly modified to satisfy gauge invariance with 
respect to some unphyical degrees of freedom introduced 
in the mapping. These modifications and the resulting 
spectrum in the D-Mott phase are described in the sub- 
sequent subsection. 

Within the GN model, the excitations are of course 
organized into SO{N) multiplets, but are further con- 
strained for the case of interest, TV = 8, by triality. In this 
subsection, we discuss the lowest-lying states, their multi- 
plet structures and quantum numbers, and give some use- 
ful physical pictures to aid in understanding their prop- 
erties. 



1. Results from integrability 

The lowest-lying excitations are organized into three 
S0{%) vector multiplets, which are degenerate due to 
triality, for a total of 24 particles. Four of the 28 global 
SO{S) generators may be chosen diagonal (to form the 
Cartan subalgebra). We will label the particles by the 
values of the four associated charges, denoted by the or- 
dered quadruplet {Ni,N2, N3,N4^), and defined by 



(4.32) 



(no sum on a). In this notation, one 5*0(8) multiplet con- 
tains the states (known as fundamental fermions) with 
only one of the four Na — ±1, and all others equal to zero. 
The remaining 16 degenerate states have Na = ±1/2 for 
a — 1,2,3,4, which are divided into those with an even 
number of Na = +1/2 (the even kinks) and the remainder 
with an odd number of Na = +1/2 (the odd kinks). The 
reasons for this terminology will become apparent later 
in this section. Each particle has a mass m and disperses 
(due to Lorentz invariance) as ei{q) — \J m? + q^, with 
momentum q. Since the electron band operators Cpia are 
defined relative to their Fermi momenta kpi, the actual 
momenta of each particle is offset from the GN model 
momentum, q, by some amount. We will return to these 
"base" momenta later in this subsection, as well as to the 
other physical quantum numbers of the excitations. 

At somewhat higher energies there is another multi- 
plet of 28 "particles" , which transform as an antisym- 
metric second rank S0{8) tensor. This multiplet can 
be viewed as two-particle bound states of the funda- 
mental Gross-Neveu fermions, or equivalently under tri- 
ality as bound even-even or odd-odd kinks. Indeed, 
of these 28 states, 24 have two zero charges and two 
Na = ±1. The other four are bound states of a fun- 
damental fermion with it's anti-particle (an excition in 
the semiconductor picture, below), so they do not carry 
any of the four quantum numbers. Each of the 28 "par- 
ticle" states has a mass 7712 = a/Sto. Finally, for energies 
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e > £ci<l) = ^\/rn?- + (?^/4, a two-particle continuum of 
(unbound) scattering states exists. 




FIG. 3: The mean field picture of the SO{N) GN model. 
There are four flavors of relativistic massive fermions V'a, 
with dispersion ei((;) = ±a/to^ + . The negative en- 
ergy bands are filled, while the positve energy bands are 
empty. As in a semi-conductor, the positive and negative 
energy bands are separated by a finite gap 2m. 



understood from the marginality of the interactions un- 
der the RG scaling transformation. The picture is that of 
a "semi-conductor", as indicated schematically in Fig. 3. 
These massive Dirac particles and their antiparticles may 
be identified with the fundamental fermion 5*0(8) vector 
multiplet. The even and odd kinks likewise arise from ap- 
plying the same decoupling to the even and odd fermion 
representations of the Hamiltonian. 

While Eq. 4.33| is correct for 5*0(00), it requires cor- 
rections otherwise. For finite iV, the chiral "order pa- 
rameter" A fluctuates around its vacuum value, and 
these fluctuations generate attractive interactions be- 
tween the GN fermions. The attractive interactions lead 
to the formation of two-particle bound states, whose mass 
ma = 2m[l - {■k'^ /2N'^) + OiN""^)] approaches twice the 
fermion mass for iV ^ 1, due to the weakness of the 
interactions in this limit. For 5*0(8), however, the inter- 
fermion interactions are not weak, and the bound states 
have the strongly reduced mass m2 = V3m. The 28 
bound states of two fermions form the abovementioned 
rank 2 tensor multiplet. A priori, one might expect three 
such multiplets, arising from bound states of the three 
sets of fermions. We will see, however, in the next section 
that this does not lead to any new particle content. In- 
deed, the particles in the tensor multiplet can be equally 
well viewed as bound states of fundamental, even, or odd 
fermions. 



2. Mean-field picture 



It is instructive to see how these excitations arise in a 
mean-field treatment of the GN interaction. The mean- 
field treatment becomes exact for the SO{N) generaliza- 
tion of the GN model for large even N . To carry it out, 
we e mp loy the Dirac fermion version of the Hamiltonian, 
Eq. 4.4. In the mean- field approximation, the bilinear 
i^X^^'i'a acquires an expectation value, and the "quasi- 
particle" Hamiltonian density becomes 



(4.33) 



where A = 2(7(-0jT^Va) is a mean- field gap parame- 
ter. The mean- field Hamiltonian is simply that of four 
massive Dirac equations. It is easily diagonized in mo- 
mentum space, using il^a{q) — exp(zri(g)T^/2)-0a, where 
ri(g) — coi''^{vq/ is), which gives 



H 



MF 



dq 
2^ 



ei(g)V'iT^V'a, 



(4.34) 



with €i{q) = -y/ w? -\- q^ and the mass m = A. From the 
diagonalized form it is straightforward to determine the 
MF estimate, 



(4.35) 



for the general SO{N) case, where A ~ i is a momen- 
tum cut-off. The exponential dependence on g can be 



3. Semi-classical picture 



These excitations can be readily understood in the 
semi-classical limit of the Bosonized Hamiltonian. In 
this language, particles correspond to classical sohtons, in 
which the phase fields 9a connect different vacuum (clas- 
sical minimum energy) values at a; = ±00. The winding 
numbers of these solitons have a direct connection to the 
5*0(8) charges, since by bosonization 



AOa = 0{oo) - 0(-oo) 
= 2n Na. 



dxdxOa 



(4.36) 



Thus each of the GN particles labelled by the four quan- 
tum numbers (A^i, N2, N3, N4) can be readily transcribed 
into a semi-classical soliton. The fundamental fermions 
are those configurations in which one of the four phase 
fields 9a changes by ±27r. The second type of soli- 
ton changes all four 9a fields by ±7r, which changes 
cos 9a ^ — cos 9a , but Icavcs the vacuum energy un- 
changed. The 2^ = 16 possible "kinks" form the semi- 
classical analog of the even and odd kink 5*0(8) octets. 
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FIG. 4: The 24 mass m excitations of the S'0(8) GN 
model, projected into the {Ni,N2) — {Q/2,Sz) plane. 
Full and open circles indicate the "fundamental" fermions 
and kinks, respectively. All 24 excitations lie on the unit 
sphere in the full four-dimensional Cartan space. The 
equivalence of a fundamental fermion and an even and 
odd kink can be seen graphically by simple vectorial ad- 
dition, e.g. the odd kink (—1, 1, 1, l)/2 and the even kink 
-(l,l,l,l)/2 add to form the GN fermion (-1,0,0,0). 



E. Consequences for the D-Mott phase 

We are now in a position to discuss the nature of 
the ground state and excitation spectrum in the D-Mott 
phase, using the technology of the GN model. 



1. Gauge redundancy 

To proceed, we must first describe the boundary con- 
ditions and gauge-redundancy needed to fully specify the 
model. Since the phase fields were originally introduced 
to bosonize the (physical) electron operators, the chiral 
electron phases are defined only moduli 2tt, 



(t)pia{x, t) ^ (t>Pia{x, t) + ^TlApiaix, t) , 



(4.37) 



where the Apia are integers. These integers describe a 
sort of gauge redundancy in the description: semiclassi- 
cal phase configurations which differ only by a different 
choice of Apia are to be treated as identical quantum 
states. Furthermore, as for any gauge theory, physical 
operators must be gauge-invariant, i.e. unchanged un- 
der the operation in Eq. 4.37, which can be performed 
locally. 



While the even and odd kinks exist for general N, for 
the special case of SO{8), triality implies that the kinks 
and fundamental solitons are on an equal footing. This 
is most easily seen using a simple graphical construction. 
Construct an A^/2-dimensional space (for N even) with 
axes 9a- In this space, draw a lattice consisting of a point 
for each vacuum configuration of the fields. All possible 
solitons can be represented on this lattice as lines con- 
necting different points to the origin (see Fig. 4). The 
fundamental fermions are then the line segments to the 
neighboring points along the axes. For N = 8, how- 
ever, there are another 16 points equidistant to the origin, 
which represent the even and odd kinks (for N > 8, these 
are further from the origin, while for N < 8 they are in 
fact closer). In this case, the even or odd kink segments 
form equally good orthonormal axes in this space, and 
viewed in this basis, the three sets of particles cyclically 
exchange their roles. 

One can also compose these particles by vectorial addi- 
tion. For instance, an even and an odd kink can be added 
to form a fundamental fermion. The 2-particle bound 
states may also be visualized in this way, by adding e.g. 
two different fundamental fermion vectors. From this 
construction it is easy to see that any such 2-particle 
state can be equally well composed from two even or two 
odd kinks. There is thus only a single 28-fold tensor mul- 
tiplet of 2-particle bound states. 



2. Uniqueness of the ground state 

From both the standpoint of the fermionic GN Hamil- 
tonian and its bosonized sine-Gordon form, the system 
appears to exhibit broken symmetry. The conventional 
GN model has a spontaneously broken "chiral" symme- 
try: the hamiltonian is invariant under the chiral trans- 
formation ■0a — !■ T^Va, however, the chiral order param- 
eter A — 2g(-!/)Jr^?/'a) and changes sign under this 
transformation. In the bosonization language, this trans- 
formation corresponds to 9a —>■ 9a + tt (for all a) , which 
takes cos 6*0 — > — cos 6*0. The bosonic model appears to 
have even more broken symmetries, i.e. there is a count- 
ably infinite set of semiclassical vacua, related by the 
additional transformations 9a 9a + ^nua, for integer 
na- 

On physical grounds, however, we expect that the D- 
Mott phase has no broken symmetry, and consequently 
a unique ground state. To reconcile this apparent dis- 
crepancy, we must account for the fact that the phases 
9a are not gauge invariant. Indeed apparently different 
vacua may represent the same physical state in a dif- 
ferent gauge. To establish the physical equivalence be- 
tween different vacua is a rather tedious and technical 
exercise, although straightforward. In Appendix ^ we 
carry through this exercise and demonstrate that all of 
the semi-calssical vacua do indeed correspond to the same 
physical state. Thus, as expected, there are no broken 
symmetries in the D-Mott phase and the ground state is 
unique. 
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3. Quantum numbers 

To connect the GN results with the physical two-leg 
ladder system, we now consider the quantum numbers 
of the various excitations. Each quantum number corre- 
sponds to some conserved quantity in the system. The 
most physically interesting are the charge, spin, momen- 
tum along the x direction, and parity (or equivalently, 
momentum in the y direction). The charge and spin are 
conserved quantities corresponding to continuous global 
symmetries, so we can work directly with the Hermitian 
generators 



Q = Jdx 4 



pja pja ' 



(4.38) 
(4.39) 



Since the translational and leg-interchange symmetries 
are discrete, we should really speak of the unitary oper- 
ators themselves. Since right- and left-moving particles 
in band j carry quasi- momentum ±kpj, respectively, the 
translation operator is simply 



(4.40) 



c 



is the total number 



where Npj = J dx'^^c 

of electrons in band j with chirality P. The quasi- 
momentum is defined by = exp(iPa;). Because 
the anti-bonding (band 1) operators have ky = tt, the 
parity or translation operator in the y direction is 



(4.41) 



where Ni = Njh + Nr,i. 

In the weak-coupling limit, U <C t,t±, there are two 
additional conserved quantities, the band spin differ- 
ence 3^2 = Sf — S2 and the relative band chirality 
P12 = {N^i - - N^2 + N^2)/2. 

It is useful to rewrite these expressions in terms of the 
bosonized phase variables. Because the symmetry gener- 
ators involve spatial integrals of fermionic bilinears, they 
can be expressed in terms of the winding numbers AS^ 
and their conjugates 

/oo 
dxdxfaix) = (Pa{oo) - fa{-Oo)- (4.42) 
-oc 



Using Eqs. 2.18-2.21, we find 



Q = 


A01 

TT 


(4.43) 




A92 

277 ' 


(4.44) 


''12 — 


A^3 

2tt ' 


(4.45) 


P12 = 


A04, 

~2^' 


(4.46) 



= exp 
Ty = exp 



Aipi 



[k 



Fl 



AOa 



- {AOi + Aifi) 



(4.47) 
(4.48) 



As discussed in the previous section, the winding num- 
bers of the da are just the 5*0(8) conserved charges. Thus 
the first four conserved quantities can be directly tran- 
scribed for all the GN excitations, i.e. {Q, S^, Sf2, P12) = 
{2Ni,N2, N3, N4). These are tabulated in the first three 
columns in Table 1. 

The momentum and parity of the particles are more 
complicated, however, because Eqs. 4.47 - 4.4^ contain the 
conjugate fields Ai^i, Aip^. As such, P^ and Py are not 
simply determined from the 5*0(8) charges Na- The ad- 
ditional physics required is the operator content of the 
original electron problem. 

To sec how this comes in, let us imagine a local opera- 
tor {{Na};x), which creates the particle with charges 
{Na} when acting on the ground state, i.e. 



O\{Nahx)\0)^\{Na};x) 



(4.49) 



where |{-/Va}) is the quantum state with one excited {Na} 
particle localized at x. Now consider the exponential of 
a phase field ipa- It can be rewritten as the line integral 
of the momentum conjugate to 0a, i.e. 



(4.50) 



where Haix) = dx^a/ {^t^) is the momentum conjugate 
to 9a, i.e. [9a{x),Ilb{x')] ~ i6abS{x — x'). Since the mo- 
mentum Ha generates translations of the phase 9a, the 
exponential operator creates a soliton of size 27rA^a lo- 
cated at the point x. In order to have the correct winding 
numbers, the desired quantum operator must thus have 
the form 

Ot({iVj;x) = A[{0J;{A^J]e-^^»^»/2. (4.51) 

Here we have included an arbitrary function A of the 9a 
fields, which does not wind the phase and thus does not 
affect the 5*0(8) charges. 

To determine A, we next impose gauge invariance. 
Consider first the operators which create the fundamental 
fermions, with only one non-zero Na = ±1. The creation 
operator takes the form 



0± = Aa[{0a}]e 



THVa+ea)/2 



(4.52) 



where we have removed a factor e^**°/^ from A to de- 
fine AJ. This is desirable because the last factor (up 
to a Klein factor Ka) is simply the GN fermion operator 
"^Jpa ii^pa for the lower sign) . Now O J must b e inva riant 
under all possible gauge transformations, Eq. 4.37. It is 
a straightforward exercise to show that the most general 
form for AJ is the same for all the fundamental fermions, 
and is given by 
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(4.53) 



where ^ ' indicates a sum over all quadruplets of integers 
with even, and 



O = gi(^i+«2+e3-e4) 



(4.54) 



Note that A does not include a term proportional to 
unity, which implies that a single GN fermion is by it- 
self not gauge-invariant and hence unphysical. Instead, 
physical particles have an attached operator Os (or its 
counterparts with extra factors from the term). Og 
represents a Jordan- Wigner "string" , and can be rewrit- 
ten only non- locally in terms of the fermion fields. It 
modifies the momentum and statistics of the fundamen- 
tal fermions to those of the physical excitations. 

It is now straightforward to determine the quasi- 
momentum and parity of the fun damental fermions using 
the translation operators in Eqs. 4.47 4.48. In particular, 
we must have 



e 



a ' 



a 



(4.55) 
(4.56) 



The left-hand sides of Eqs. f4.55| - p^ can be evaluated by 
employing the commutators of the Bose fields to obtain: 



TxOaTj. 
Tx<faT~^ 



kF2)Sa 



(4.57) 
(4.58) 



and 



Tydafy^ ^Oa+ 2^(5„4, (4.59) 
fy^aT-^ = ^a + 2nSal. (4.60) 

Eq. 4.55 can be written as a product of three terms: 

f^O^f-^ = (E') ^ TxOst-' 

xj;^eT^(¥>.+ea)/2j,^-i^ (4.61) 

Consider the first term. Using the above commutators 



IS mvari- 



one can readily show that the sum in Eq. 4.53 
ant under x— translations. The string, however, carries 
momentum ~ tt: 



TxOsT- 



(4.62) 



This momentum must be added to the "bare" momentum 



of the GN fermion, given by the last term in Eq. 4.61 
A similar calculation for Ty shows that is again in- 
variant, but Os carries transverse momentum tt. The 
resulting net momenta of the fundamental solitons are 
summarized in the last two columns of Table. 1. 

Similar manipulations hold for the kink excitations. In 
particular, the even kink creation operators must obey 



O. 



e± 



{ka} 



X e 



and similarly for the odd kinks. 



o± 



{fea} 



,Ti(v°-e°)/2 



(4.63) 



(4.64) 



Note that the choice to fa ctor out t he right-moving even 
and GN fermions in E qs. 4.52 , 4.63 and left-moving odd 
fermions in Eq. 4.64 is arbitrary. A right-mover can 



always be converted to a left-mover and vice-cersa by 
absorption of an e^'^ factor into a redefinition of the 
"string" part of the soliton creation operator. For the 
even and odd kinks, the above choice is particularly 
convenient, since the right-moving even fields and left- 
moving odd fields are exactly bare electron operators 
and hence manifestly physical. We see from Eqs. |4.63- 



4.64 that the kinks have the quantum numbers of bare 
electrons. A remarkable consequence of this result is that 
the number of single-electron degrees of freedom has ef- 
fectively doubled relative to the free fermi gas, since each 
of the 16 kinks can be created with arbitrary momentum 
(relative to its base momentum of ±kpi), including par- 
ticles beiow the former fermi sea and holes above it. The 
momenta calculated from Eqs. 4.63| - 4.64 complete the 
last two columns of Table 1 . 



iNi,N2,N3,Ni) 


Q Px Py 


(1,0,0,0) 
(0,1,0,0) 
(0,0,1,0) 
(0,0,0,1) 


2 

1 TT TT 
TT TT 
2kpl 


(l,l,l,l)/2 
(l,-l,-l,l)/2 
(l,l,-l,-l)/2 
(1,-1,1, -l)/2 


1 1 kpl TT 
1 — ^ kpl TT 
I \ kp2 Q 
1 -i kp2 


(l,l,l,-l)/2 
(l,-l,-l,-l)/2 

(l,l,-l,l)/2 
(l,-l,l,l)/2 


1 ^ —kpl TT 
1 — ^ —kpl TT 

1 i -kp2 

1 -i -fc^2 



Table 1: Physical quantum numbers of the mass m par- 
ticles labelled by their four U{1) charges. The antiparti- 
cles are obtained by changing the sign of all the quantum 
numbers. 



F. SU{2) invariance and spin multiplets 

We conclude this section with a remark on SU (2) in- 
variance and the excitations with spin. In Table 1, we 
have classified the mass m excitations in the D-Mott 
phase by the four U{1) charges {SO{8) Cartan gener- 
ators) Na- Another almost equivalent choice is to la- 
bel the particles by their charge, momentum, total spin 
= s{s + 1) and spin projection . It is fairly trivial 
to relabel the kinks in this way, reflecting their natural 
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correspondence with single-particle spin-1/2 excitations. 
They may be grouped into doublets with s = 1/2 and 
S"^ = ±1/2, e.g. the (l,l,l,l)/2 and (1, -1, -1, l)/2 
kinks form a spin-1/2 doublet with Q = le and P = 
(kpijir). For the GN fcrmions, this is less trivial. The 
(±1,0,0,0) solitons are spin zero, and correspond to 
charge ±2e Cooper pairs with zero momentum. Simi- 
larly, the (0, 0,0,±1) fermions carry neither charge nor 
spin, and may be regarded as dressed particle-hole pairs 
carrying only momentum. The remaining four solitons 
are more nontrivial, however. Their spin content may be 
brought out by re-fermionizing the total spin operator, 
as in Eq. [4.17, 



S = j dx [J r{2 



with chiral currents. 



(4.65) 



(4.66) 



with A,B,C — 3,4,5. Thus the three Majorana fields 
77p^, with A = 3, 4, 5, transform as a triplet of spin s = 1 
operators, i.e. 



— tt l/pc- 



(4.67) 



The other 5 GN Majorana fields commute with S, and 
hence represent spin-singlet operators. 

As was shown in the previous subsection, the physi- 
cal GN excitations consist not of GN fermions but rather 
required an attached string Og- Fortunately, the string 
does not carry any spin, i.e. 



O.SOl = S. 



(4.68) 



Thus the true soliton excitations (GN fcrmions+strings) 
satisfy the same transformation rules with respect to spin 
as the bare Majorana fermions. The four remaining soli- 
tons (0, ±1,0,0) and (0,0, ±1,0), which involve the four 
Majorana Fermions r]^ with A = 3,4,5,6, can therefore 
be decomposed into an s = 1 triplet and a spin-zero 
singlet (both with Q = and P = (7r,7r)). The triplet 
can be regarded as a minimum energy magnon, while the 
singlet can be grouped with the (0,0, 0,±1) solitons as 
another particle-hole excitation. 

With the SU{2) invariance realized, we can tabulate 
the particles in the GN model in a slightly different way, 
classifying them by their SU{2) multiplet (i.e. s = 0, 
1/2 or 1), charge, and momentum. To label the parti- 
cles classified in this way, the abelian {Ni, N2, N^, N4) 
notation is no longer convenient, since it does not re- 
spect the SU (2) invariance. Instead, we can schemati- 
cally indicate the 8 particles in the vector multiplet by 
iJa and the 28 in the tensor multiplet by tjaVb (remem- 
bering that rfsr/A — — t]^'?!?)- For convenience, we list the 
8 GN fermions and the 28 mass \/3m bound states in this 
way in Table 2 (we do not list the remaining 16 particles, 
since they have the quantum numbers of electrons and 
are easily remembered in this way). 
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1 


±{kpi ~ kp2) 


TT 



Table 2: Physical quantum numbers of the mass m 
(above horizontal line) and mass y/im (below horizon- 
tal line) particles. 



CORRELATION FUNCTIONS AND 
PHYSICAL PROPERTIES 



We have seen that in the weak-coupling limit, the two- 
leg ladder possesses an enhanced symmetry. The effec- 
tive theory in this limit is the 50(8) GN model, which 
is both exactly integrable and exhibits a remarkable "tri- 
ality" . In this section we will discuss a variety of the 
resulting physical consequences. 

The most remarkable consequence of triality iSp- the 
equality of the single-particle and two-particle gapsEja 
The 16 kinks have the same quantum numbers as the bare 
electrons at the former Fermi surface. The single-particle 
gap, defined as the minimum energy needed to add an 
electron or hole to the system, is thus simply Ai = m. 
The 8 GN fermions, however, have charge Q = ±2 or 
Q = Q, and thus represent excitations corresponding to 
an even number of electrons and/or holes. For instance, 
electron or hole pairs can be added with zero net momen- 
tum via the (±1, 0, 0, 0) solitons, while spin-1 excitations 
may be added with momentum (tt, tt) via the (0, ±1, 0, 0) 
and (0,0, ±1,0) solitons (more precisely the 775 state). 
The gap for all 8 minimal energy 2-particle excitations is 
also A2 = m. 

The equality of the single-particle and two-particle 
gaps is in marked contrast to the behavior of other more 
conventional insulators. In a band insulator (such as the 
two- leg ladder at half- filling with t±^ ^ t), the single- 
particle gap is just the band gap, while the two-particle 
gaps are twice as large: A2 = 2Ai. Another famil- 
iar case is the strong-interaction limit U ^ t. In this 
case, the single particle gap is huge Ai ~ U, while the 
lowest two-particle (e.g. spin) gaps are much smaller 
A2 ~ t'^/U ^ Ai or indeed vanishing (A2 = 0) for or- 
dered or quasi-long-range ordered antiferromagnets (e.g. 
d = 2 or single-chain Hubbard models). The detailed 
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mathematical mapping between the GN, odd, and even 
fermion fields allows us to extend the relationship be- 
tween the single- particle and two-particle properties be- 
yond the values of the gaps, as we detail below. 

First, we will discuss several correlation functions 
which characterize the spin and charge dynamics of the 
system. The most interesting of these are the single- 
particle spectral function, measurable by electron photo- 
emission, and the dynamic spin structure factor, which 
is probed by inelastic neutron scattering. Other inter- 
esting correlators include the current-current correlation 
function, which determines the conductivity, the density- 
density correlation function and pairing correlation func- 
tion, which can be measured in numerical simulations. 



A. Single-particle spectral function 



First consider the single-particle Greens function. 



Gl(fc,T)=^, 



''T^r(0|ai+,Jx,r)aL(0,0)|0), 



(5.1) 



where is the (Euclidean) time-ordering symbol and 
fc = (kxyky). The right-hand side is independent of a 
by spin-rotational invariance; however, we choose to de- 
fine the spectral function for fixed a, i.e. no sum is im- 
plied above. In general, the single-particle spectral func- 
tions can be extracted from Gi by Fourier transforma- 
tion. Defining Gi{k^iLu) — J dr Gi{k,T) exp^iur), one 
finds 



1 



ImGi(fc, iuj ^ uj + iS) ^ip(fe, oj) + Aih{—k, —uj) 



(5.2) 

where the particle and hole spectral functions are 
Aip(fc,c^) = ^|(n|aljO)|'j(fc-fc„)'5(cc>-i;„), (5.3) 

71 

A,hik,u) = Y,\H'^ia\0)f - K)Siu; - E^). (5.4) 

n 

Here we have abbreviated S{k) = 2TT5{kx)Sky.o- The 
task is then to evaluate Gi(fc,r). In the weak-coupling 
limit studied here, this is obtained in two stages. We 
first integrate out the electron fields cl^{k) , c^^{k) , for 
I A: — kpi\ > A, which can be accomplished perturbatively 
in the interactions, since the energy denominators are 
finite away from the Fermi momenta. The perturba- 
tive corrections to the free-electron Gi(fc,T) are there- 
fore small in these regions. Within the cut-off region of 
width 2A, we must employ the full RG treatment. The 
RG scales the problem onto the GN model, which thus 
applies at the lowest energies. 



For the electron spectral function, the non-interacting 
result ^ip//j is trivial, since single electron states are ex- 
act eigenexcitations. The result is 



A'ij,{k,uj) = S{uj-e,{k))9{uj), 



A%{k,u;)^S{iu + ei{k))9{uj) 



(5.5) 
(5.6) 



where ei(k) = —tcoskx — -j-cosfcy. Interactions of 
course modify this form somewhat, leading to some spec- 
tral weight away from lo = ei(k), and a broadening of 
the delta-function peak in Aip/h for some momenta. In 
weak-coupling, away from small lu, however, both effects 
are small. We will return to them after we consider the 
behavior of the spectral function for small frequencies. 

The low-frequency limit of Aip is dominated by mo- 
mentum near the Fermi points. Transforming to the 
slowly- varying Luttingcr fields, we have 

G,iPk,, + q,ky,-T) « / da;e-'«-(c,,„(a;,r)4„(0,0)). 



(5.7) 

for g ^ 1, where kyi — {2 — i)n. Unfortunately, integra- 
bility does not give exact forms for the time-dependent 
correlation functions in Eq. A considerable amount 
can be learned, however, from the exact excitation spec- 
trum and from approximate methods. 

The spectrum determines the support of j4ip/^(k, w). 
From Eq. 5.3, it is clear that ^]^p//j(k, cj) is non-zero 
only when there exists an excitation (or more than one) 
with momentum k, spin s = 1/2, charge ±e and en- 
ergy UJ. From Table 1, we see that the sixteen kinks 
have exactly the appropriate quantum numbers for all 
possible momenta near the four (i.e. two pairs) fermi 
points with either — ±1/2 and charge ±e. We there- 
fore expect that Aip^ii(k,uj) first becomes non-zero for 

1^. Since the kinks 



Lu = ^y{kx - kpiY + {ky 
are isolated particles with a fixed energy-momentum re- 
lation, these excitations give a sharp delta- function peak 
in Ap/fi. It is natural to identify this peak as the contin- 



uation of the non-interacting delta-function in Eqs. 5.5- 
5.6 to the region near the Fermi points. At higher ener- 
gies other states should contribute to the spectral weight. 
A quick consideration of the quantum numbers is suf- 
ficient to conclude that none of the mass ^/Sm bound 
states have the appropriate quantum numbers (e.g. all 
have charge zero or ±2e). The next lowest- lying excita- 
tions with the quantum numbers of individual electrons 
are in fact scattering (unbound) states of a kink and a 
GN fermion. For instance, a (1,1,1,— 1)/2 kink and a 
(0,0,0,1) fermion can form a scattering state with the 
quantum numbers of a spin up electron with momentum 
{kpi + (7,7r), with q 1. Similarly, other combinations 
of kinks with the (0, 0,0,±1) GN fermions contribute 
to the single-particle spectral weight at (±3kpi + (?, tt), 
{±kp2 + q, 0), and {±3kf.2 + Q, 0). All these form continua 
with Lu > edq) — l^m^ -\- {(l/2)'^ at each momenta, since 
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the energy at a particular momenta can always be in- 
creased continuously by shifting the kink and the bound- 
state momenta in opposite directions. Further excitation 
of more than one (0, 0, 0,±1) quanta leads to spectral 
weight at all momenta separated by an even multiple of 
2kpi, i.e. k = {{2n + l)kpi,Tr). The excitation gap for 
such a point increases, however, by an additional factor 
of m as each GN ferniion (i.e. factor of {2kpi,0) away 
from the Fermi points) is added. Furthermore, the higher 
harmonic contributions to the spectral function are ex- 
pected to have small amplitudes, as they involve multiple 
scatterings of the original injected electron (see below). 




-|-m2). (5.9) 




FIG. 5: Schematic plot of the single-particle electron 
spectral function Aip(kx,ky — n,uj). The curve below 
the continuum indicates a sharp resonance, i.e. a delta- 
function peak in Aip, which accquires a finite width once 
it passes inside the continuum due to the onset of decay 
processes. The continuum above energy 2m coincides 
with the creation of (unbound) scattering states of a sin- 
gle kink and a GN fcrmion. 

To understand the magnitude of yl]^p//j(k, w) in the al- 
lowed regi ons requires a knowledge of the matri x ele ments 
in Eq. 5.3, or of the full Green's function in Eq. 5.7. Since 



exact results are unavailable for these quantities, we con- 
sider instead the mean-field approximation. Without loss 
of generality, let us consider momenta near the particu- 
lar Ferm i point k » (fcj-i, tt) and spin = -1-1/2. Using 
Eq. 5.7 and Eqs. 4.28 - 4.31 , the bare electron operators 
can be rewritten exactly in terms of even Fermion fields, 

(c«iJx,r)4iJ0,0)) = {rr,,{x,T)rn\{0,0))- (5.8) 

In the mean-field approximation, the exact eigenstates 
are created by the rotated operators for particles ■0Ha(fc) 
and holes iplai^) (^^^ IV. D. 2), so that the above ex- 
pectation value can be computed by a simple rotation. 
One finds 

</(fc,i+g,^,^) = 



The g-dependent factor out front arises from the rotation 
to the ip spinor, and is an alogo us to a "coherence factor" 
in superconductivity. Eq. |5.9| captures a simple and ap- 
pealing physical effect. Although single-kink states ex- 
ists for all momenta, their contribution to the spectral 
weight is suppressed for g <C — m by the "coherence fac- 
tor" above. Such a negative momentum corresponds to 
the addition of an electron at a momentum inside the 
Fermi sea. Interactions deplete the Fermi sea slightly 
near the Fermi points, allowing electrons to be added, 
but with a weight that vanishes as g ^ — oo, i.e. deep 
within the sea. Similarly, the hole spectral function has 
weight outside the Fermi sea, since some excited parti- 
cles exist as part of the ground state. Unfortunately, the 
mean-field approximation is not sophisticated enough to 
capture the continuum ai uj > Cc{q), since the ip fields cre- 
ate exact eigensta tes h i this limit. Thus the continuum is 
not present in Eq. 5^. On physical grounds, however, we 
expect that it will be similarly suppressed for momenta 
inside the former Fermi sea. 

Having determined the behavior of Aip/f^ near the 
Fermi points and momenta separated from them by 
even multiples of (2kpi,0), we return to the question 
of the spectral weight away from these points, at en- 
ergies away from the non-interacting single-particle en- 
ergy ei(k). In the non-interacting system, states with 
one added electron plus additional neutral electron-hole 
pairs in fact form a continuum away from the single- 
particle energy. Consider, for instance, ky = tt. The 
lowest energy continuum of states with a single added 
particle-hole pair (plus one electron) consists of those 
states in which both electrons are infinitesimally above 
the Fermi surface (fci = k2 = kpi) and the hole makes 
up the missing momentum (fca = 2kpi — k). This begins 
at e2{kx, 7r) = cos(2fcj^i — kx) — t±/2, which is below the 
single-particle energy (i.e. £2 < £1) for kpi < kx < 3kpi, 
crossing zero again at kx — 3kpi. It does not contribute 
to the single-particle spectral function, however, due to 
orthogonality. In an interacting Fermi liquid, we would 
expect that an added electron can scatter into these 
states (i.e. emit a low energy particle-hole pair), and 
some weight would appear in Aip^j^ associated with the 
continuum states. 

In the weakly-interacting ladder, an added electron 
away from the Fermi points can also scatter to create neu- 
tral excitations, and some weight will appear in the re- 
gions near the non-interacting continuum for uj ^ e2(k). 
This continuum away from the Fermi points should merge 
smoothly into the continuum above lu — 2m at the 
Fermi points (and the higher harmonics, e.g. {±3kpi, tt)). 
Glearly, since the single-particle energy begins below the 
continuum near the Fermi momenta and it is above the 
continuum far away, it must cross into the continuum 
at some point. Where it is below the continuum, we 
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expect that the spectral function retains a sharp delta- 
function peak. Once it passes above, however, the single- 
particle mode can decay into the continuum states, and 
should acquire a small width. Putting this behavior to- 
gether with the spectral function near the Fermi points 
and higher harmonics, we arrive at the schematic single- 
particle spectral function illustrated in Fig. 5. The most 
dramatic feature is the sharp delta-function peak near 
the Fermi points which crosses into the continuum and 
acquires a width at higher energies. 



B. Spin structure factor 

The spin spectral function can be defined in a similar 
way to the single-particle one. Consider the structure 
function 



U'x •' 

xTr{0\Slix,T)SliO,0)\0), 



(5.10) 



where the lattice spin operator is S£{x) = '^la'^^ep- 
The spin spectral function is obtained from this in 
the usual way 

-Im S^^ {k,iuj ^ Lj + iS) As{k,uj)S^\ for a; > 0. 

TT 

(5.11) 

The spectral decomposition of As is 
^(k, ^) = J E I (0|^>) r '^(k - K)S{lo - E,,). (5.12) 

n 

As for the single-particle case, we expect that the spin 
spectral function will be approximately equal to its non- 
interacting value for uj m. A straightforward calcula- 
tion for the non-interacting problem gives 



a,b—±l 2 




h cos \/ 1 — r'^ + r sin + at± /2 

b cos — y 1 — r'' — r sm ha — cos fc 

2 2 2 



(5.13) 



where r = + a^{cosky - l))/(2sin^). Clearly, 
the non-interacting spin spectral function is considerably 
more complex than its single-particle counterpart. This 
is because the neutral spin-one excitations in the non- 
interacting system are particle-hole pairs, and thus com- 
prise a continuous spectrum. The result in Eq. 5.13 is 
plotted in Fig. 6. 
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FIG. 6: Intensity plot of the non-interacting spin spectral 
function at (a) ky = 0, (b) ky — tt. The darkness is pro- 
portional to the spectral weight, white indicating regions 
of phase space in which no particle-hole pairs exist. 



For ky 



0, the particle and hole must come from 



the same band. In this case low energy excitations exist 
near kx = 0, when both are taken near the same Fermi 
point, and near kx = ±2fcFi {kx = ±2kp2 are equivalent 
to these values modulo 2n), when the particle and hole 
are taken from opposite Fermi points (still in the same 
band). Two branches of continua (from the bonding and 
anti-bonding bands) form as the momenta of the parti- 
cles and holes are varied. For ky ~ tt, the particle and 
hole must come from opposite bands. In this case, the 
kx — states have an energy of exactly t± , since this re- 
quires a vertical transition. Low energy excitations exist 
near kx = ±{kp2 — kpi), due to particle-hole pairs taken 
from two right-moving or two left-moving Fermi points. 
They also exist near kx = tt, when the particle and hole 
are taken from opposite sides of the Fermi surface. Ex- 
tending these two branches of excitations gives the form 
shown in Fig. 6(b). 
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FIG. 7: Intensity plot of the interacting spin spectral 
function at (a) ky = 0, (b) ky = ir in the presence of 
interactions. In the low-energy portion, various gaps de- 
velop and excitations with sharp delta- function peaks are 
present (see Sec. V.B). Note that the minimum energy 
spin excitations occur at k = (tt, tt). 

As for the single-particle spectral function, introduc- 
ing interactions allows for additional structure, and the 
low-lying excitations are raised up to energies of order m. 
In particular, from Table. 2, we see that the lowest lying 
neutral triplet states are the mass m GN fermions ?y3,4,5 at 
k = (tt, tt). The next highest energy neutral triplets are 
the mass bound states. The ?73,4,5??7,8 have momenta 
{±{kFi - A;F2),7r), while the ?73.4,5f?4,5.3 and ry3,4.5??6 have 
momentum (0, 0). Above these exist continua dispersing 
like ec{q) — + away from all the aforemen- 

tioned points and {±2kpi,0) (the excitations at these last 
points arise from certain pairs of unbound kinks) . Since 
there are no sharp resonances (delta-function peaks) in 
the non-interacting limit, we expect that the mass m and 
^/3m peaks must broaden at higher energies to merge into 
the continua found there. A schematic form is shown in 



Fig. 7. 



C. Optical conductivity 

Another quantity of considerable experimental rele- 
vance is the optical conductivity. We are interested in 
the real part of the conductivity, defined by 



Re a{uj) = Im 



Il{iuj Lo + iS) 



(5.14) 



where the (k = 0) current-current correlator is 

n{iLu) 1 dxdre''^'' {TrJi{x,T)Ji, {0,0)) . (5.15) 

The electrical current operator is 

Je{x) = ^ (^al{x)af,{x + I) ~ al{x + l)a^(a;)^ . (5.16) 



To evaluate Eq. 5.15| , only the slowly varying (k = 0) 
component of the current is needed. Decomposing 
the latti ce f ields into their continuum components us- 
ing Eq. ^.2| and then applying the bosonization and 
refermionization rules, one finds the long-wavelength 
form 



Ji ~ sin kpidx 



sin kpi 



27r 



(G--G-). 



(5.17) 
(5.18) 

(5.19) 



From this form, the current operator clearly excites 
the mass ^/im 771772 bound states, as well as higher en- 
ergy continuum scattering states with energies above 2m. 
Since Eq. ^.14 is nothing but 1 /lo times the spectral func- 
tion of J, the zero temperature optical conductivity is 
thus zero for < ^/Zm, has a sharp (delta-function) peak 
at u) — ^/Sm, and a threshold with continuous weight for 
uj > 2m. Based on the mean-field picture, we expect the 
spectral weight in the two-particle continuum to have a 
square root singularity (due to the van Hove singularity 
at the bottom of the band - see, e.g. Ref. Es), i.e. 



a{uj) w A5{uj - \/im) + B 



UJ - 



2m 



e(w- 277i). (5.20) 



At non-zero temperatures it is diff icult to determine 
<j{uj) from the Kubo formula, Eq. 5.14 . Instead, the gen- 
eral features can be argued on more conventional trans- 
port grounds, borrojaiing heavily from recent results of 
Damle and Sachdevcj for spin dynamics of gapped 2- 
leg ladders. The important physical effect for T > is 
the presence of a non-zero equilibrium concentration (oc 
g-m/T-j q£ activated excitations. In the semi-conductor 
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analogy, these are activated particles and holes. In prin- 
ciple, all 24 mass m states have identical equilbrium con- 
centrations; for T ^ m we expect that we can neglect the 
much smaller activated densities (0(e~^™/"^)) of bound 
states. 

In this case the low-frequency conductivity can be es- 
timated using a simple Drude argument. We focus on 
the charged species of mass m, i.e. the 771,2 fundamen- 
tal fermions and the kinks. Each of these contributes in 
parallel to the conductivity a term of the Drude form. 



Co 



1 + iujT 



(5.21) 



with (To = nr/m, t a scattering time, and n the density 
of thermally excited carriers. Using the Boltzmann dis- 
tribution, we have n ^ ^/mTe~"^^'^ for T <^ m. The 
average separation between particles 1{T) = 1/n is thus 
much larger than their typical wavelength A(r) ^ 2ii /p ^ 
l/\/mT, obtained by equipartition. The particles thus 
behave essentially classically except during a collision, 
when they scatter strongly (as known, e.g. from the 
exact S-matrices for the SO{N) GN model), and their 
scattering time is determined simply by the time be- 
tween collisions: t ~ l{T)m/p ~ T^^e"^/"^. The ex- 
ponential dependences in the dc conductivity thus can- 
cel, and (To ^ 1/^/mT (the same result is obtained from 
the Einstein relation erg = f^I-^)- principle, the di- 
mensionless numerical prefactor in this relation could be 
obtained using the methods of Rcf. |4^, but we content 
ourselves here simply with the scaling form. Note that 
although the height of the Drude peak diverges as T ^ 0, 
its width shrinks much more rapidly (exponentially) , and 
the weight at w = is negligible at low temperatures. 

Turning to the higher-frequency features (for oj ~ -s/Sm 
and CO 2m), we expect that scattering between the in- 
jected bound states or particle-hole pairs and the ther- 
mally excited carriers will occur on the same character- 
istic timescale, r. These peaks therefore also acquire ex- 
ponentially small widths (0(1/t)) for T ^ A. Since the 
overall spectral weight in Re 17(0;) must be conserved, we 
expect the heights of these features to diverge much more 
strongly than the Drude peak, i.e. a{uj — ^/Sm, T) ^ t ^ 



and ct(w = 2to, T) - ~ e 



m/2T 



Indeed, the op- 



tical conductivity presumably satisfies universal scaling 
forms near these points, i.e. 



Re {(7(0))} 
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FIG. 8: Optical conductivity at finite temperature. At 
low temperatures, all the features become exponentially 
sharp, with a width Slo ^ e"'"/'^. In this limit, the "exci- 
ton" peak at w = y/3m retains a non-zero weight, sharp- 
ening into a delta-function, and the peak near oj = 2m 
is also exponentially high. By contrast, the Drude peak 
at w = has vanishing weight at low temperatures, its 
height diverging only as erg ^ 1/^/mT. 



D. Equal-time spatial correlators 

Numerous other correlators can be measured at equal 
times in numerical simulations, and sometimes experi- 
mentally (e.g. static structure factors). The properties 
of essentially any such correlator can be deduced from 
the GN spectrum, as summarized in Table 2. Due to the 
Lorentz invariance of the GN model, intermediate states 
with a finite energy e give rise to exponentially decay- 
ing spatial correlation functions with the corresponding 
length = v/e. 

For completeness, we quote two examples here. The 
pair-field correlator, defined by A{x) — ai^{x)a2i{x) has 
the correlation function 



(A(x)A^(O)) ~ Aie^™l"l/"- 



+^36 



■ A4e 



—i2kFix 



(5.23) 



^ ' \y/TY.^[{uj -2m)T\ |w-2m|<TO ' 

(5.22) 

for T <^ m, where and Ej are universal scaling func- 
tions. A schematic illustration of the optical conductivity 
at finite temperature is given in Fig. 8. 



for \x\ ^ 1. Here the first term comes from the mass 
m "Cooper pairs" and the second from the correspond- 
ing bound states. In the prefactors to the exponentials, 
Aj, we have neglected sub-dominant spatial dependences, 
which generally will have power law forms. For exam- 
ple, due to the one-dimensional van-Hove singularity for 
adding a pair above the threshhold energy m, one ex- 
pects Ai{x) ~ for large x. Similarly, the real-space 
density-density correlation function is 
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B3 + B4{-lY cos(fcj.i - kp2)x 



'V3m\x\/v _|_ . . . (5.24) 



Here ni{x) = alai^)'^eai^) ~ 1- The real-space spin-spin 
correlation function has an identical form, except with 
B2 = 0. 



VI. GENERIC INTERACTIONS AND SO(5) 
SYMMETRY 

In the previous sections, we focused on the properties 
of the D-Mott phase which occurs with generic predom- 
inantly repulsive interactions in the two-leg ladder. In 
weak coupling this phase exhibits a remarkable 5*0(8) 
symmetry with dramatic physical consequences for both 
two-particle and single-particle properties. As remarked 
in the introduction, there exists an SO(h) subalgebra of 
the full 50(8) group whose vector representation "uni- 
fies" superconductivity and antiferromagnetism. Thus all 
the consequences of this 5*0(5) symmetry, proposed by 
Zhang as a phenomenological model for the cuprates, are 
shared by the D-Mott phase. A number of authors have 
proposed exactly 5'0(5)-invariant lattice models includ- 
ing, in a recent paper by Scalapino, Zhang and Hanke 
(SZH), a two- leg ladder model.CJ SZH derived a complex 
phase diagram for this model in the strong coupling limit 
in a space including both repulsive and attractive inter- 
actions. In this section, we will develop the necessary 
technology and study in weak-coupling both the SZH 
model aud-ather generic 50(5)-invariant two-leg ladder 
systems EJ'EJ 

In fact, the focus on 50(5)-invariant models is less re- 
strictive than might be naively expected. Indeed, in our 
numerical studies of the full RG equations at half-filling, 
all weak coupling two-leg ladder models we have studied 
{including those with attractive interactions ) scale under 
the RG onto the 50(5)-invariant manifold. Within this 
manifold, we have observed five attractors, including the 
D-Mott phase and a Luttinger liquid (C2S2) phase con- 
tinuously connected to the non-interacting Fermi liquid, 
in which all the elementary excitations remain gapless. 
The remaining three attractors are all massive phases, 
and comprise an S-Mott phase similar to the D-Mott 
phase but with approximate "s-wave" pair correlations, 
a ciarge-density-wave (CDW) phase with a density-wave 
at k = (T^jTr), and a spin-Peierls (SP) phase without a 
density wave but with kinetic energy modulated at the 
same wavevector. We group the D-Mott with the latter 
three to form four dominant phases. We will see that (in 
weak coupling) while all of these dominant phases share 
Zhang's 50(5) symmetry, each one possesses a distinct 
higher 50(8) symmetry. The different 50(8) symme- 
tries are related in rather simple ways which have ramifi- 
cations for the critical points between the different dom- 
inant phases. 



A. SO{5) Symmetry 



We begin by reviewing some basic properties of the 
5*0(5) symmetry, demonstrating in the process the re- 
lation to the 50(8) symmetry already discussed. The 
50(5) algebra was originally designed to rotate the five- 
component vector containing the real and imaginary 
parts of the D-wave pair field and the three components 
of the staggered magnetization. A set of operators which 
performs this function was introduced by ZhangciJ - these 
are the 10 generators of 50(5), which are conveniently 
grouped into the antisymmetric matrix 





Qp 

ReH^r 
ReH^ -ImH^ 





-ImH^ 
-ImH,, 







~Sy Sx 



(6.1) 



where A,B = 1 ... 5 spans the matrix of generators and 
Kab — —Kba- The various components are defined as 
bi-linears in electron operators, 

Qp = ^EH(^K(fc)-l). (6.2) 

k 

S=\Y.ai{k)cT^pap{k), (6.3) 

k 

n= i^0fca„(-fe + Ar)(CTyCr)„/5a^(fe), (6.4) 



where TV = (tt, tt) is the "nesting" vector and Qp is the 
charge measured in the number of pairs of electrons rel- 
ative to half-filling. Here aa(k) is the Fourier transform 
of the lattice electron operator and a spin sum is im- 
plicit. To see that the 50(5) symmetry is just the first 
subgroup of 50(8), we need simply take the continuum 
limit of Eq. |6.l| . Straightforward but lengthy decomposi- 
tion of the lattice electron fields into their slowly varying 
components (as in Section II) and consequent bosoniza- 
tion and refermionization (as in Section IV) gives the 
exceedingly simple result 



J dx (or + Gi 



(6.5) 



where the Gp^ are precisely the 50(8) generators in- 
troduced in Section IV. Since only the five-dimensional 
upper- left block of the full matrix of SO (8) generators 
enter in Eq. |6.5| , the 5*0(5) symmetry rotates the first 
five Major ana fermions 77^, A = 1 ... 5. As discussed in 
Section IV, the first five components of this vector rep- 
resentation contain both the pair field {'^ ^71, 2) and the 
staggered magnetization ('^ 773,4,5). 
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B. Microscopically 50 (5)-invariant models and 50(5) 
spinors 

We now turn to a discussion of microscopically SO{5)- 
invariant ladder Hamiltonians. A particular example 
is the SZH model, which is the most general 50(5)- 
symmetric two-leg ladder Hamiltonian with nearest- 
neighbor hopping and only intra-rung two-body interac- 
tions. The interaction terms on each rung of the ladder 
take the form 

H^nt = U^{{ne^ - ^){nei - i)} 

e 

+V{ni-l){n2-l) + JSi-S2, (6.6) 

where £ = 1,2 refer to the two legs. S0{5) symme- 
try requires a single constraint on the three couplings: 
J = 4(f7 + V). For t/, y » t,tj_ the hopping t can 
be treated perturbatively, and SZH have determined the 
(quite complex) phase diagram in the U — V plane. 
With the weak-coupling RG, we can attempt to com- 
plete the phase diagram by exploring the opposite limit, 
U,V ^ t,t±. (One can hope to determine the behav- 
ior at intermediate coupling U,V ^ t, t± by interpola- 
tion.) Further, we can explore the generic behavior of 
other weakly-interacting 50(5)-invariant two-leg ladder 
systems which contain, for example, inter-rung interac- 
tions. 




FIG. 9: The folded Brillouin zone for the S0{5) spinor. 
The allowed momenta are chosen to be in the grey area. 
For a two-leg model, the transverse momentum ky takes 
two values 0,7r. In our convention, only ky = excita- 
tions are inside the folded Brillouin zone. 

To do so, we need a means of constructing 50(5)- 
invariant models in weak coupling. For lattice roadels, 
such constructions have been discussed by HenleyEj and 
Rabello et. al.,LJ and apphed by SZH to the two-leg 
ladder. Here we generalize these methods to the chiral 



fermions operators which appear in the linearized con- 
tinuum model obtained in the weakly-interacting limit. 
Since the Hamiltonian is built from electron operators, 
we need to introduce spinor representations of 50(5). Wa 
begin with the lattice construction of Rabello et. al.,E2l 
defining the four-component spinor as 



aa{k) 



(6.7) 



where AT = (tt, tt) is the nesting vector of the Fermi sur- 
face. To avoid double-counting, the allowed momentum 
k in the spinor only runs in the "folded" Brillouin zone, 
whose size is half of the original one, as shown in Fig. 9. 
For the two-leg ladder model, in which the only trans- 
verse momenta are ky = 0, tt, it is possible to specify 
the folded Brillouin zone by simply setting fcj, = in the 
above spinor. The factor 0^ (which in the general two- 
dimensional case is a non-trivial function with absolute 
value one) can be taken to be unity with this convention. 
Re-expressing the spinor in terms of the band electron 
operators and Fourier transforming gives 



*(x) = /dfc,*(fc„0)e*'=^" 



(-l)"cL(x) 



(6.8) 



In the continuum limit valid for weak coupling at low 
energies, a chiral decomposition is possible: 



with chiral spinors defined by 



(6.9) 



(6.10) 



was can- 



To obtain Eq. |6.10| , the (-1)=" factor in Eq. | 
celled using the relation kpi + kp2 = tt. 

The advantage of the spinor basis over the electron 
band operators Cpia is that they transform simply under 
50(5). In particular, under a unitary transformation 
generated by the operator 



eyip{i9ABKAB), 



3.11) 



where A, B ~ I . . .5^ the spinors ^'p transform according 
to 



3.12) 



where the spinor indices a,b = 1 ... 4, and T{9) = 
exp(i6'^Br^^) is the rotational matrix for a spinor. Here 
^ i[r^,r^]/4 where the F-* are five generalized (4 
by 4) Dirac matrices, discussed in detail in Appendix |y. 
They satisfy the usual Clifford algebra 



{r^,r-} = 25^ 



(6.13) 



Using the spinors, we can break down all fermion bilin- 
ears into irreducible representations of 50(5), i.e. gener- 
alized currents. Three "normal" sets, which involve one 
and one \E' spinor, carry net momentum zero or (tt, tt): 
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7„ = ^ 

'-' P — ^ pa ^ pa^ 



^ ,T,t 



(6.14) 
(6.15) 
(6.16) 



The three currents in Eqs. 6.14-3.16 transform as an 
5*0(5) scalar, vector, and rank 2 antisymmetric tensor, 
respectively. A second set of currents (and their her- 
mitian conjugates) appear "anomalous" , and carry net 
momentum (±2^:^2,0) or (t: ±2kp2,T^) '■ 



I^P = ^' p,{RT-)ab^ Pb- 



(6.17) 
(6.18) 



These two currents, which transform as a scalar and a 
vector under 50(5), require the introduction of the ma- 
trix 



R = 



1 

-1 



(6.19) 



where 1 is the two by two identity matrix. Note that it is 
straightforward to show that the matrices i?r^^ are sym- 
metric, so that a non- vanishing anomalous tensor current 
cannot be defined. A simple counting verifies that the 
above set of currents completely spans the space of elec- 
tron bilinears. There are 1-1-5-1-10 = 16 currents in 
Eqs. 3.14-6.16, and an additional 2 • (1 + 5) = 12 currents 



in Eq. 3.17-6.18 and their complex conjugates, for a total 
of 28 = 8 • 7/2 independent bilinears. 

In weak coupling, we must generically consider all Her- 
mitian products of two bilinears which are (1) invariant 
under 50(5) and (2) conserve quasi-momentum. Ne- 
glecting purely chiral terms (which, as in Sec. 13, only 
renormalize velocities), there are then five allowed cou- 
plings. The interaction Hamiltonian density takes the 
form 



'Hint — 9sJrJl + QvJrJl + QtJa^ Jl 

h, i IrIl + h-c. '>+hJ Itif + h.c. 



(6.20) 



Note that momentum conservation forbids forming a 
quartic interaction from one normal and one anomalous 
current. 

The above Hamiltonian represents the most general 
50(5) invariant ladder theory with weak interactions. 
The five coupling contants {gs, gv, gt, hs, hy) specify a 
five-dimensional manifold within the more general nine- 
dimensional space of U{1) x SU{2) symmetric theories. 
This manifold is determined explicitly by a set of linear 
equations, given in App. ^, which relate the five 50(5) 
invariant couplings to the 9 t/(l) x SU{2) couplings which 
were introduced in Sec. II. Because the 50(5) manifold 
posesses higher symmetry, it closes under an RG trans- 
formation. The five RG equations describing the flows 
within the 50(5) manifold are given explicitly in App. ^. 



C. SZH Model and Four Dominant Phases 

The weak coupling phase diagram for the SZH model 
can now be obtained, by numerical integration of the 
50(5) invariant RG flow equations (Eqs. |El[-p^). The 
initial (bare) values of the 50(5) coupling constants 
are obtained by taking the continuum limit of the SZH 
model. For each initial set of bare parameters, the phase 
is determined by bosonizing those couplings which grow 
large under the RG transformation, as described in Sec. 
II. The resulting weak coupling phase diagram is shown 
in Fig. 10. 

Four new phases appear in addition to the D-Mott 
phase which occurs for predominantly repulsive interac- 
tions. In the "C2S2" region in Fig. 10, all five cou- 
plings scale to zero under the RG. This Luttinger liquid 
phase thus retains all the gapless modes (2 charge and 2 
spin, hence C2S2) of the original non-interacting electron 
system, and thereby has (an approximate) chiral 50(8) 
symmetry. 
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FIG. 10: Phase diagram of the 50(5) symmetric SZH 
model plotted in the U-V plane with J = 4{U + V) and 
U,V<^t = t^. 

We group the other three states together with the D- 
Mott as four dominant phases. In the S-Mott phase, 
the interactions diverge in the same way as in the D- 
Mott case, given by Eqs. 3.2, with the modification 
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In the 50(5)-invariant 



notation, this corresponds to changing the sign of hs 
and hy. Semiclassically, the only change in the behav- 
ior is that {(pp-)s-Mott = (Vp-)D-Mott + TT = TT. The 
6a± and 9p+ fields are unaffected, so that the S-Mott 
phase still has short-range pairing. It is, however, of ap- 
proximate s-wave symmetry, with AiA , > due to the 
shift in ipp-, as can be seen from Eq. 3.8. It is inter- 
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esting that the strong-coupHng "s-wave" paired state on 
a rung, | tij— ) + | — iTi) is identical in the ladder leg 
and band bases, and corresponds to an on-site pairing or 
singlet state. In contrast, pairing in the strong-coupling 
D-Mott state is across the rung of the ladder, as depicted 
schematically in Fig. 11. 

In the SP and CDW phases, the ratios of diverging 
couplings are somewhat different. In particular, fo^^i^ii 
are irrelevant and 



4/r2 = (T)6?2 = (±)J&?2 



'12 



(±)2<i = .g > 0, 



(6.21) 
(6.22) 



where the upper and lower signs hold in the SP and CDW 
phases, respectively. These modifications imply a fairly 
dramatic change in the behavior relative to the D-Mott 
and S-Mott states. In fact, the SP and CDW are dual 
to the D-Mott and S-Mott, respectively, in the follow- 
ing sense: each is obtained from its dual counterpart 
by interchanging (pa-- and 9a-- ■ Because of this inter- 
change, the pair fields fluctuate wildly even locally, and 
CDW = 0. Instead, these two phases 
break discrete Z2 symmetries. 

To explore this in detail, consider the order parameters 



Be{x) = ^[aljx + l)af„(x) + h.c] 
riiix) = aljx)aea{x) - 1. 



(6.23) 
(6.24) 



The field B({x) is the local kinetic energy, while ne{x) 
is the local electron density relative to half filling. The 
two order parameters differ in symmetry since Bi{x) is 
even and ni(x) is odd under a Z2 particle-hole symme- 
trya^^(x) (— l)^+'^a|^(a;). Using the usual relations to 
rewrite Be and ne in terms of chiral operators, bosonizing, 
and applying the semi-classical results (common to both 
the SP and CDW phases) (6ip+) = (6i<,+) = (ipa-) = 0, 
one obtains 



1 



(Bfix)) ^ (-l)-+^(cos(-^,_)), 



(neix)) ^ (-l)-+^(sin(-^,_)) 



(6.25) 
(6.26) 



Since {ipp-) = Q,tt in the SP and CDW phases, re- 
spectively, we find {Be)sp ^ (— 1)"^+^, {ne)sp — and 
{Be)cDw = 0, (n^)cDW ~ (-1)"+^. The SP phase thus 
breaks only the discrete symmetry under translations by 
one lattice spacing (the translation in the y-direction is, 
of course, equivalent to parity), while the CDW phase 
breaks both translational symmetry and the Z2 particle- 
hole symmetry. These broken symmetries can be de- 
picted easily in the strong coupling limit, as shown in 
Fig. 11. 
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FIG. 11: Schematic illustration of the four dominant 
phases, drawn for simplicity in the strong-coupling limit. 
In the D-Mott and S-Mott phases, neighboring rungs 
contain essentially decoupled pairs. Adjacent rungs are 
highly correlated in the SP and CDW phases, which fur- 
thermore break parity symmetry. 



D. SO(8) Symmetries and Degeneracies of the 
S-Mott, SP, and CDW Phases 



Since the four dominant phases appear on essentially 
equal footing, one might suspect that the S-Mott, SP, 
and CDW phases exhibit 5*0(8) symmetries similar to 
that of the D-Mott phase. We shall see that this is in- 
deed the case, but that the 50(8) algebras are different 
in each state. 

Consider first the S-Mott. In the previous subsection, 
it was shown that the S-Mott is related to the D-Mott 
by a TT shift in (pp^ . It follows that if we define 



a = 1,2,3 
a = 4 



(6.27) 



and (p^ = fai the bosonized Hamil tonian in the S-Mott 
phase takes the form of Eqs. ^^^^ (with Oa, <<5a replaced 
by 0^, (/j^). Consequently, a refermionization into the GN 
form is again possible. In particular, upon defining 



'Ipa I p^^ 



A^l...Q 
A = 7,8 



(6.28) 



the S -Mo tt Hamiltonian takes the "canonical" GN form 
(Eq. |4.6 ) in terms of the ?7^^. The sign changes in 
Eq. |6.2g imply that the generators of the 50(8) sym- 
metry in the S-Mott phase (are different from those of 
the D-Mott phase. For instance, = — whose 
spatial integral is not an 5*0(8) generator in the D-Mott 
case. However, since the Majorana fermions for the two 
phases are equal for A = 1 ... 6, the D-Mott and S-Mott 
do share a common 5*0(6) subalgebra. 

Similar constructions can be performed for the SP and 
CDW phases. Recalling the duality of the previous sub- 
section, we choose 
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Oa a = 1,2 

tpa- a — 3 
ipp- a = 4 



fa a =1,2,4 

0rr- a = 3 



(6.29) 



for the SP phase. Similarly for the CDW, we take 



qCDW 




a = 1,2 
a = 3 
TT a = 4 



(6.30) 



and (^a^^ = (^5^^. As before, the GN form can be 
retained. The appropriate Majorana fermions in these 
cases are 



SP 

'IPA 



j Vpa 



A = 1...5,7,i 

A = 6 



and 



CDW 

'IPA 



rjpA 
PriP 



A 
A 



1...5 

6,7,8 



(6.31) 



(6.32) 



Like the D-Mott and S-Mott, the SP and CDW phases 
share a common 5*0(6) symmetry. Moreover, the D-Mott 
and SP share an SO{7) symmetry, as do the S-Mott and 
CDW. 

A final calculation is possible now that the appropriate 
bosonizcd variables have been established. In Sec. [V and 
App. the uniqueness of the D-Mott ground state was 
established. We also expect a unique ground state for 
the S-Mott phase, but have yet to establish it. In the SP 
and CDW phases, discrete symmetries are broken, and 
one expects at least a two-fold degeneracy in the ther- 
modynamic limit. Using the techniques applied earlier 
(gauge equivalence of semiclassical solutions) to the D- 
Mott, we can determine these degeneracies. Details can 
be found in App. ^ The result of such an analysis is 
that the S-Mott indeed has a unique ground state, while 
the SP and CDW ground states are each exactly two-fold 
degenerate. 



E. Full Set of SO{5) Fixed Points 

Because of the relative simplicity of the 5*0(5)- 
invariant manifold (5 coupling constants versus 9 for the 
general case), it is possible to perform an exhaustive de- 
termination of the possible asymptotic scaling trajecto- 
ries under thcRG. To do so, we insert the power -law 



ansatz of Eq. |2.26| into the 50(5) RG equations, Eqs. ^A 
E5. This set of five coupled algebraic equations can 
be solved exactly, in contrast to the corresponding set 
of nine U{1) x SU{2) equations, which have proved in- 
tractable. One finds fourteen solutions, delineated in Ta- 
ble 3. Five represent the states discussed so far: the gap- 
less C2S2 and four dominant 5'0(8)-symmetric phases. 
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Table 3: Fourteen algebraic solutions of the 50(5) RG 
equations 

Of the remainder, five represent critical points. Con- 
sider first the D-Mott^S-Mott transition. Taking the 
values in Table 3, one finds that semi-classically the fields 
6'p-|_ , , and 6^^ are pinned as in the D-Mott and S- 
Mott states, but that neither the 6p- nor the Lpp- field 
appears in the interaction Hamiltonian. There is thus 
a single gapless (central charge c — 1) bosonic mode 
at the critical point. That this is indeed the critical 
point between the D-Mott and S-Mott phases can be 
seen by perturbing slightly away from the scaling tra- 
jectory. If the perturbations are small, those terms in- 
volving the gapped degrees of freedom will have negligible 
effect, and we need only include the couplings involving 
the p— fields. As argued in Sec. 0, cos nOp- terms are 
not allowed by translational invariance. The low-energy 
Hamiltonian density near the critical point (after inte- 
grating out the massive fields) thus has the form 



'^D^s = [{dxfp- f + {dcc9p-f] -Xcosifip-. (6.33) 



For A > 0, the semiclassical minimum occurs for ipp^ = 0, 
describing the D-Mott phase, while for A < 0, the mini- 
mum shifts to (pp- — TT, yielding the S-Mott phase. We 
expect that the general form of this low-energy critical 
Hamiltonian will remain valid even in strong coupling, 
though the Luttinger stiffness and velocity of the critical 
(fip- mode will shift in this case. What are the critical 
properties of this transition? The correlation length ex- 
ponent is determined by the scaling dimension of cos (pp^ . 
In a general strong-coupling situation, this is a contin- 
uously variable exponent. In weak-coupling, however, 
it is determined. In particular, refermionization implies 
cos(pp- ^ V4''"^V'4i which acts as a Dirac mass and has 
scaling dimension one. In this limit then, the correla- 
tion length ^ ^ l-^l""; with — 1. Both in strong and 
weak coupling, the dynamical exponent z = 1, as deter- 
mined by the quadratic bosonic kinetic energy. This type 
of c = 1 continuously-variable critical point is known as 
a Gaussian model, as shown in Fig. 11. Of course, in 
neglecting the massive modes, we have thrown out ad- 



25 



ditional universal physics in the weak-coupling limit. In 
particular, these massive modes have a large S0{6) sym- 
metry, which can be seen by rewriting the critical inter- 
actio n Hamiltonian density using Table 3 and Eqs. 



Dll 



A.B=1 



(6.34) 



The full weak-coupling critical symmetry is thus U (1)_r x 
U{1)l X 5*0(6). It may seem surprising that this critical 
point has lower symmetry than the massive dominant 
phases, which enjoy S'0(8)-invariance. This is a result 
unique to the weak-coupling limit. With stronger interac- 
tions, corrections to the weak-coupling scaling will break 
the S0{8) symmetry, while leaving the C/(1)_r x U{1)l 
critical symmetry (which results from truly infinite- 
wavelength physics) intact. 

Having understood the D-Mott<->S-Mott transition, it 
is clear that the SP<-^CDW transition is essentially iden- 
tical. The Hamiltonian in this case differs only via the 
interchange of 9a-- and ipa-- , which in any case are mas- 
sive at this critical point. 

The next two critical points are somewhat different. 
For concreteness, consider the D-Mott^SP transition. 
As before, three of the bosonic fields are massive, in this 
case 9p+,9cr+,9p-. These can be integrated out, leav- 
ing the fj— fields critical. However, here hoth 9^- and 
the dual ipa~ appear, so a semi-classical analysis is not 
tenable. Instead, we refermionize this single remaining 
bosonic field and its interactions after integrating out 
the three massive bosons (i.e. setting them to their semi- 
classical minima). The reduced Hamiltonian density in 
this case is 



n 



int.,GfT. 

D^SP 



(6.35) 



Here g is the finite coupling along the scaled RG tra- 
jectory, and A is a deviation from the critical trajectory 
similar to A in the Gaussian model above. Since g is 
non-zero, the rjp^ Majorana fermion acquires a gap, and 
only the single rjpa Majorana fermion is gapless at the 
critical point. This is a central charge c = 1/2 critical 
point, which uniquely identifies it as an Ising transition. 
Indeed, the Ising nature of this transition is very phys- 
ical, given the discrete Z2 parity symmetry broken in 
the SP phase. This also explains the duality between the 
D-Mott and SP phases found earlier: this duality is noth- 
ing but the usual Kramers- Wannier duality of the Ising 
model. As before, we expect the Ising critical behavior 
to be robust to corrections to the weak-coupling RG, so 
these transitions should be in the same universality class 
even with strong interactions. In the weak-coupling limit, 
the massive degrees of freedom again have higher sym- 
metry, in this case including the rjp^ Majorana fermion 
coming from the a— fields. The full weak-coupling criti- 
cal theory is thus Z2 x SO{7), where the Z2 theory is the 
conformally-invariant Ising model, as indicated in Fig. 
11. 



Not surprisingly, the S-Mott^CDW transition is also 
of the Ising variety. The "multi-critical point" in Table 3 
describes the case when all four phases come together at 
a point, i.e. when two transition lines cross. It is sim- 
ply a direct product of the two critical theories above, 
i.e. a Gaussian model and an Ising theory. It is pos- 
sible that these theories actually become coupled if one 
re-introduces interactions that were irrelevant at the non- 
interacting Fermi fixed point, but we do not explore this 
possibility here. 

The remaining four "fixed points" of the 5*0(5)- 
invariant RG describe more exotic situations. We have 
not observed any microscopic Hamiltonians attracted to 
these phases, but some of these may perhaps occur for 
some choices of bare interactions. We suspect that these 
"phases" are unstable, and hence spurious for physically 
relevant situations. Nevertheless, we discuss them briefly 
for completeness. They are mo st easily u nderstood by us- 
ing the representations in Eqs. D10| - D14 . The form of the 
S'0(5) X 50(3) case is then easily seen from the interac- 
tion Hamiltonian density (taking the positive sign for hs 
for simplicity) 
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SO(5)xSO(3)GN ^ 9 
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A,B = 1 



GAB ryA 
7^ ^ T 



GAB ryA 
7^ ^ T 



(6.36) 



These interactions are precisely those of an 50(5) x 
50(3) GN model. Both of the constituent GN mod- 
els are massive, so this represents another gapped phase. 
The solution with the opposite sign for can be con- 
verted into the same form by the canonical transforma- 
tion ?7„e ^Vrb, so it is also a gapped phase of this 
sort. The remaining two phases can be understood simi- 
larly. Note that in these cases only the scalar interactions 
gs and are non-zero. This implies that the first five 
and last three Majoranas are decoupled. Furthermore, 
in this case since gv — gt — 0, the first five Majorana 
fermions are non-interacting. They comprise a gapless 
5*0(5) Wess-Zumino-Witten (WZW) model with central 
charge c — 5/2, while the last three Majorana fermions 
combine to form an 5*0(3) GN model. The interaction 
Hamiltonian density is thus simply 



njint. 

'^SO{5)WZWxSO{3)GN 
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GAB r^A 
I, r 



(6.37) 



VII. DOPING THE D-MOTT PHASE 

In this Section we briefly consider the effects of dop- 
ing away from half-filling in the two-leg ladder. In the 
weak coupling limit of interest and with nearest neigh- 
bor hopping in the kinetic energy, the ladder at half- 
filling was argued to scale onto the soluble Gross-Neveu 
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model which posseses an exact global 5*0(8) symmetry. 
Generally, doping away from half-filling will break down 
this large symmetry, leaving only charge and spin conser- 
vation, with the much smaller U{1) x SU{2) symmetry. 
This can already be anticipated for the non-interacting 
problem: When the Fermi energy moves away from zero 
(half-filling) the Fermi velocities in the bonding and anti- 
bonding bands will in general become unequal due to 
curvature in the energy /wavevector dispersion. For weak 
doping, however, this effect is small. Indeed, in the rela- 
tivistic model derived in Section II where the dispersion 
was linearized about the Fermi points, the small curva- 
ture is ignored entirely. In the following we focus on this 
very low doping limit (x = 1 — ?i <C 1), where the dif- 
ference between the two Fermi velocities can be safely 
ignored. We thus continue to employ the linearized rela- 
tivistic model. Nevertheless, as we shall see, even within 
this limit doping away from half-filling breaks down the 
global 50(8) symmetry of the Gross- Neveu model, al- 
though in a rather straightforward manner. 

To dope we consider adding a chemical potential term 
to the Gross-Neveu Hamiltonian, H, with i/^ = H — ^Q, 
where Q is the total electron charge. This charge can be 
written 

Q = 2jdx i^l,,p„, + V'liV'.i) = 2 (Gf + Gf). 

(7.1) 

Since Q is a global 50(8) generator, it commutes with 
the full interacting Hamiltonian: [Q, H] = 0. Thus, even 
for /i 7^ the states can still be labelled by Q, which, 
along with all the generators with A,B = 3 ... 8, 
remains a good quantum number. The 50(8) multiplets 
will of course be split by the presence of a non-zero /i, 
lowering the energy of positively charged excitations and 
raising the negatively charged ones. 

The splitting of the 50(8) multiplets can be conve- 
niently visualized in the large-N "semiconductor" pic- 
ture. Of the four fermionic particle/hole excitations 
of the fundamental GN fermions, only the first one is 
charged and is shifted in energy. Specifically, employ- 
ing the semiclassical notation, the excitations (±1, 0, 0, 0) 
carry charge Q = ±2, and are shifted by an energy 
AEf^ = =f2m, as depicted schematically in Fig. 12. Pro- 
vided this shift is smaller than the energy gap, —m < 
2/1 < m, the ground state remains unaltered: The nega- 
tive energy "valence" bands remain filled and the "con- 
duction" bands empty. In terms of electrons, the ladder 
remains at half-filling. 

Similarly, the energies of the 16 kink excitations are 
split into 8 with energy m + n and 8 with energy m — /z. 
Of the 28 (two-fermion bound) states with energy y/3m, 
16 are neutral and unshifted by the chemical potential. 
Of the others, 6 have charge 2 and are shifted up in en- 
ergy by 2/i and the other 6 down by —2fi. 

The situation is more interesting when 2/i > m. In 
this case, the energy of the states in the "conduction" 



band for Cooper pairs drops below zero, and the ground 
state will be radically altered. In the large-N limit the 
new ground state will consist of filling up the negative en- 
ergy states with a Fermi-sea of Cooper pairs, as depicted 
in Fig. 12. For N — 8 the pair excitations will not be 
describable in terms of free Fermions, but one still antic- 
ipates the general picture of a conducting sea of Cooper 
pairs to remain valid. Since the Fermions interact for 
finite N, this conducting sea will be more correctly de- 
scribed as a Luttinger liquid. In the limit of very low 
doping, however, the Cooper-pairs will be very far apart 
and well described in terms of hard-core Bosons or free- 
Fermions. In this limit, the Luttinger liquid parameters 
should approach those of Free fermions. It is probable 
that the N=8 Gross-Neveu model remains integrable even 
in the doped case, since the states can still be labelled 
by the same good quantum numbers, so that exact state- 
ments about the doped Mott-insulator can be made. In 
the following we are less ambitious, using known results 
from integrability for the undoped case to infer the be- 
havior in the very low doping limit. 




y^a- Vp- 

FIG. 12: Mean-field picture of the doped 50(8) GN 
model. Only the first band (with Q = 2) is shifted by the 
chemical potential. For 2^ > m, Cooper pairs are added 
to the original ground state. 



A. Excitations with One pair present 

When 2/1 > TO the energy is lowered by adding Cooper 
pairs to the system. Here we consider first the case 
2/1 = TO + 0+, so that the concentration of pairs, denoted 
x/2, is infinitesimal. In this limit it is sufficient to con- 
sider the properties in the presence of a single Cooper 
pair. The presence of even this one pair modifies the 
spectrum of other excitations - such as the spin or single- 
particle gaps - as we now briefly discuss. 

Consider flrst the spin-gap, i.e. the energy of the 
lowest-lying spin s = 1 excitation. In the undoped case, 
the lowest-lying triplet states are the r/3,4,5 fundamen- 
tal fermions, with momentum (tt, tt) and energy to. As 
known from integrability, these excitations interact via 
an attractive interaction with the other GN fermions, in- 
cluding the single Cooper pair present due to doping. In- 
deed, the binding energy is known exactly and given by 
Eb = {2 — -\/3)to. With the single Cooper pair present, 
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an s = 1 magnon can be added to the system into a 
bound state with the Cooper pair, costing a reduced en- 
ergy m — Eb = (%/3-l)m. Thus the spin-gap at infinites- 
imal doping a; = 0+ is reduced from the undoped value 
of m down to (\/3 — l)m. There are, of course, also un- 
bound s = 1 excitations which can be created well away 
from the Cooper pair, with energy m. In fact, for x — > 
the energy m spin excitations will dominate the spectral 
weight. The spectal weight for the lower energy s = 1 
bound states will presumably vanish linearly with x. It is 
worth emphasizing that the discontinuity in the spin-gap 
at infinitesimal doping a; = 0"*" is a general feature due to 
the presence of a magnon/Cooper-pair bound state in the 
undoped Mott insulator, and is not an artifact of weak 
coupling. If such a bound state survives strong coupling, 
as suggested by numerical RG on the two-leg ladder, a 
discontinuity should be present. 

It is also instructive to consider the energy gap for 
adding single electrons in the presence of the single 
Cooper pair. Adopting a convention where Q > corre- 
spponds to "hole" doping, consider the energy to add a 
single electron with charge —1. A single electron can be 
created by adding a kink excitation, for example an even 
kink with (— 1,— 1,— 1,— 1)/2 in the semiclassical nota- 
tion. When ^ = this costs an energy m, but is shifted 
up in energy for non-zero chemical potential: Ei = m-\-fi^ 
as depicted in Fig. 13. When 2/i = m + and the single 
Cooper pair is added, the energy to add the electron can 
be lowered from Ei = (3/2)to by binding the kink to the 
(1,0,0,0) Cooper pair. This forms a charge Q = 1 hole 
state: an odd kink with (1, — 1, — 1, — 1)/2. The associ- 
ated binding energy equals m, as follows directly from 
triality (at fj, = 0). Thus, at infinitesimal (hole) dop- 
ing the energy to add an electron drops by m, down to 
to/2. As for the case of the spin-excitations, one expects 
a continuum of unbound single electron excitations, at 
energies above 3to/2. 



B. Excitations with Many pairs 

For 2/z > TO the "conduction band" for Cooper pairs 
will be partially occupied. In this case, one expects a 
continuum of low energy particle/hole excitations cre- 
ated by exciting pairs across the Fermi "surface". For 
the SO{8) Gross-Neveu model the Cooper pairs will pre- 
sumably interact with one another, so that the semicon- 
ductor picure of a non-interacting Fermi sea will not be 
quite correct. Rather, the spinless gas of Cooper pairs 
will presumably form an interacting Luttinger liquid. In 
any event, one expects a continuum of low energy excita- 
tions in the Cooper pair fluid, presumably with a linear 
dispersion relation. One might hope that the velocity of 
this mode as a function of doping x might be accessible 
from integrability of the doped Gross-Neveu model. 

It would also be very interesting to study the en- 
ergy of the spin-one excitations with finite doping. A 



s = 1 magnon added to the system will interact via 
an attractive interaction with the sea of Cooper pairs. 
For infinitesimal doping [x = 0+) the corresponding 
spin-gap energy was lowered due to the formation of 
a magnon/Cooper-pair bound state. With many pairs 
present, this energy will presumably be further lowered, 
as depicted schematically in Fig. 12. 

Finally, we consider briefly the spin one excitations 
at energies above threshhold. These states would con- 
tribute to the spin-spectral function, accessible via in- 
elastic neutron scattering in the doped ladder. Gen- 
erally, we expect a continuum of states above thresh- 
hold, corresponding for example to adding a magnon 
at (tt, tt) and simulatenously exciting multiple "particle- 
hole" pairs in the (Cooper-pair) Fermi-sea. This con- 
tinuum should contribute to the spin-spectral function 
at any given momentum. For example, at momentum 
(tt, tt), multiple particle- hole pairs with zero net momen- 
tum will contribute spectral weight at all energies above 
threshhold. Due to this continuum of states, we do not 
expect any delta-functions in the energy dependence of 
the spin-spectral function in the doped ladder. 

This expectation runs contrary to arguments put for- 
ward by Zhang for the existence of a sharp tt— resonance 
in the spin-spectral function in the superconducting 
phase of models which exhibit an exact SO (5) symmetry. 
Zhang's argument has recently been applied to the doped 
(power-law) superconducting phase of the two-leg ladder 
by Scalapino, Zhang and Hanke.Cj Below, we briefly re- 
consider Zhang's argument for the sharp tt— resonance, 
and show that in addition to S0(5) symmetry, it relies 
on the existence of a condensate in the superconducting 
phase. Being one-dimensional, however, a true conden- 
sate does not exist in the "superconducting" phase of 
the two-legged ladder. In our view, this invalidates the 
argument for a sharp delta-function tt— resonance in the 
doped ladders, even in the presence of exact S0(5) sym- 
metry. 

Zhang's argument for the tt— resona nce rests on the 
fact that the tt operators, defined in Eq. 6.4, being global 
S0(5) (and S0(8)) generators, are exact eigen-operators 
even with non-zero chemical potential: 



(7.2) 



where = H — and the subscript a labels the 
three components of the tt— operators. This implies that 
for any eigenstate H^\E) = E\E) with energy E, the 
triplet of states Ii-a\E) are also exact eignestates, but 
with energy E+2fi, provided they are non- vanishing. De- 
note the exact ground state of the doped ladder with N 
Cooper-pairs as |iV), which satisfy Hf^\N) = 0. Adding 
an additional Cooper-pair is accomplished with the oper- 
ator Ot{x) = Os{x)i;l{x). The zero momentum Fourier 
transform of this operator, Oi{k — 0), creates a state 
with + 1 pairs, which can be decomposed as. 



O+{k = 0)\N)=c\N + l) + ..., 



(7.3) 
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where the dots denote excited states with N + 1 pairs 
present. Following Zhang, we can use the tt— operators 
to rotate Cooper-pairs at zero momentum into a triplet 
of magnons at momentum (tt, tt), since 



(7.4) 



with a = 3, 4, 5. Acting with the H— operator on Eq. |7.3| 
and using the above commutation relation and the fact 
that Os commutes with IIq, one obtains, 



Oaik = 0)\N)^^na\N- 



1) 



(7.5) 



where Oa{x) = Os{x)ria{x). The left hand side is a spin 
1 triplet of states with momentum (tt, tt), built by adding 
a magnon to the N-pair ground state. Due to the SO (5) 
symmetry, the states on the right side, XlajA^ + 1), are 
exact eigenstates with energy 2fi. As argued by Zhang, 
the equality between the left and right sides implies that 
the triplet of magnons will contribute a delta-peak in the 
spin-spectral function at energy 2fi - the tt— resonance. 

However, this conclusion rests on the assumption of 
a non- vanishing overlap between Oa\N) and IlaliV -I- 1). 
But in the thermodynamic limit, the squared overlap, 
|cp, is simply the (Bose) condensate density, since 

c= (iV + l|C'j^(fc==0)|A^). (7.6) 

While non-zero in a 2d superconductor, for the two-leg 
ladder the condensate density is zero, and the argument 
for a delta-function tt— resonance is invalid. The vanish- 
ing condensate density is a general property of Id systems 
which follows from the Mermin- Wagner theorem in the 
thermodynamic limit. For finite N at fixed pair density, 
we expect c to decay like an inverse power of the system 
length L. 
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FIG. 13: Spin-gap as a function of doping x. The spin- 
gap is discontinuous at a; = due to the formation of 
magnon-Cooper-pair bound states. 

For this reason, we expect that a finite length 5*0(5)- 
invariant ID model exhibits a S{uj — 2fi) peak in the 
spin spectral function at momentum (tt, tt) with weight 
(coefficient^ decreasing as some power Zhang has 

suggestedcS that the spin spectral function may have a 



corresponding algebraic singularity in frequency in an in- 
finite system. To address the fate of this finite-size peak 
in the thermodynamic limit, we consider now an approx- 
imate calculation of the doped spin-spectral function in 
the infinite system. 

To this end, we must determine the GN operator con- 
tent of the lattice spin operator. Using the techniques of 
Sec. [V, it straightforward to show that the decomposi- 
tion of S*^ (x) contains a term 



(7.7) 



Of course, many other operators are also present, but give 
either negligible or identical contributions to the spec- 
tral function in the regime of interest. In the D-Mott 
phase, the string Os has negligible qualitative effects in 
correlation functions, since the bosonic fields 6a are all 
locked (i.e. only weakly fluctuating) around 0a = in 
that case. In the doped system, however, there is an im- 
portant modification to the 9i field. Since the d erivative 
of this field is just the pair density (Eq. 4.43| ), its av- 
erage value has a mean slope {6i{x)) = —2kfx, where 
Kp — 7r(l — Ti)/2 is the Fermi wavevector for the sea of 
Cooper pairs; recall that 1 — n is the concentration of 
holes in the system. Furthermore, there will be fluctua- 
tions of 01 around this mean value, corresponding to the 
density and phase waves of the Cooper-pair fluid in the 
CISO state. 

To account for both these effects, we redefine 0i{x) 
—2kpx + Oi{x), treating the shifted (zero mean) di field 
as a free Bose field, as appropriate for a free Fermi or 
Luttinger- liquid system. The remaining three {02-, 03, 04) 
fields remain locked, and we therefore set these to zero 
inside the Jordan- Wigner string Os- This gives 



S+{x) 



(7.1 



As carried out for the undoped case in Sec. 0, the spin 
spectral function can be extracted from the analytically 
continued Fourier transform of the imaginary time spin- 
spin correlation function 



Su,{x,t) = {S-{x,t)S+ (0,0)). 



(7.9) 



Using Eq. 7^, one then finds 



xV',^(a;,r)V'L(0,0) 



■l{eiix,T)-ei{o,o)) 



(7.10) 



To proceed, we require a calculation of the above ex- 
pectation value. The simplest natural approximation, 
which will be our first attempt, is to decouple the charge 
(1) and spin (2) sectors, calculating the 0i correlator as 
appropriate for a Luttinger liquid (i.e. from a free Bose 
theory) and the ip2 correlator using the "semiconductor" 
free-fermion operators. In particular, one finds 
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-K/4 



(7.11) 



where K is the Luttinger parameter of the Cooper-pair 
fluid; K = 1 corresponds to free fermions, as is appro- 
priate for very low dopings. Here we have set the Fermi 
velocity of the Cooper pair sea to one. The fundamental 
fermion correlator is approximately 



(V;,„(a;,r)VL(0,0)) 



MF 



dp 
2^ 



e'P''-''^P'>''Q{T), (7.12) 



where Q is the Heavyside step function. To simplify 
Eq. 7.12, we have neglected to include the mean-field "co- 
herence factors". Because these are non-singular, their 
neglect only modifies the final result by an overall smooth 
momentum-depe ndent amp litude factor. Multiplying the 
two terms in Eqs. 7.11-7.12, performing the Fourier trans- 
form and analytically continuing to real frequencies gives 
the spin spectral function 



f 

Im / dx dp dr ■ 



-ipx+{Ltj~ei {p+k—KF)+i5)T 



(a;2 + r2)^/4 



-e(r), (7.13) 



where 5 = 0^ is a positive infinitesimal. Singular behav- 
ior can only arise from the large x^t power-law behavior 
of the denominator. For large a;, the oscillating exponen- 
tial implies that the integral is dominated by p « 0, so 
that the dispersion ei can be linearized around this point. 
Doing so, the p and x integrals can be readily performed. 
Up to an overall constant prefactor, one finds 

poo 

AfF(7r-yfc,7r;w)~Im / r-^/2g(^-ei(fc-K.)+.5)r^ 
Jo 

(7.14) 

This integral can be related to a Gamma function by an- 
alytic continuation. Carrying this out carefully gives the 
final mean-field result 

Af^{Tr-k,n;u;) - \u;~ei{k- K,FT^+-^e[uJ~eiik - kf)]. 

(7.15) 



As suggested above, Eq. 7.15 indeed exhibits an alge- 
braic singularity. For momentum (7r,7r), k = above, 
and the Fermi-level condition eiin) = 2/i for the Cooper- 
pair fluid indeed implies the singularity is located at 
Lv = 2/i, identifying it with the putative "Pi resonance" . 
Note, however, that within this approximation identical 
"resonances" appear at all momenta, including a lower 
energy one ed k — up- Moreover, the resonance becomes 
more singular when the Luttinger parameter K decreases 
approaching a delta function as K ^ 0, whereas the Pi- 
resonance should approach a delta function in the oppo- 
site limit oi K ^ oo where the Cooper-pair fluid devel- 
ops off-diagonal long-ranged order. Thus, it is unclear 
whether the above resonance for the Id model has any 
connection with the two-dimensional Pi-resonance. 



Moreover, further reflection on the nature of the mean- 
field approximation used above, leads us to question 
the validity of the singular behavior at finite frequency. 
While it might well be correct for the 0{N — oo) GN 
model, the fundamental fermions, e.g. Cooper pairs and 
magnons, are strongly interacting for the N = 8 case of 
interest, as evidenced, e.g. by the 0(1) binding energy 
for the mass bound states and the degeneracy of the 
fundamental fermion and kink excitations in the D-Mott 
phase. While interactions will not significantly modify 
the 6i correlator above (since the Cooper-pair fluid re- 
mains a Luttinger liquid), they would appear to have a 
drastic effect upon the '02 Green's function. In general, 
this Green's function describes the propagation of a sin- 
gle massive injected particle into and interacting with a 
Luttingpi liquid. Similar problems have been extensively 
studiedjy and one finds that the massive particle will 
generally radiate both energy and momentum into the 
Luttinger liquid, decaying in the process. From such de- 
cay processes, we generally expect a finite lifetime and 
hence broadening of the algebraic singularity above. For 
large N, the interaction and hence the broadening would 
be small, but we see no reason for this to be the case 
for N = 8. Furthermore, one might naively expect that 
the minimum energy singularity at k — Kp would sur- 
vive, since it is at the bottom of the V'2 band and thus 
naively has no states to decay to. The mean-field approx- 
imation, however, misses the existence of bound states, 
including the Cooper pair-magnon bound state which lies 
below the band minina at very low doping. In general, 
we expect that even the k — Kp fundamental fermion can 
decay into this bound state (radiating excitations in the 
Luttinger liquid) in the interacting system, washing out 
the algebraic singularity even here. 

In summary, the above argument suggests that above 
threshhold the spin-spectral function at finite doping will 
be smooth as a function of energy, with no singularities. 
Since this conclusion is based on a number of physical 
arguments and approximations, we cannot rule out the 
possibility of some high-energy singular structure. Cer- 
tainly singular behavior at a; = 2/i would be a remarkable 
phenomenon. On a much firmer standing is the spin-gap 
threshold energy, which is presumably a universal func- 
tion of doping X for the Gross-Neveu model. One might 
hope that the precise functional form for this energy gap 
is accessible via integrability. 
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APPENDIX A: RG EQUATIONS 

For the weakly interacting two-leg ladder with particle- 
hole symmetry at half-filling there are nine marginal non- 
chiral four fermion interactions, as discussed in detail 
in Section II. The leading order renormalization group 
(RG) flow equations for the corresponding nine interac- 
tion strengths are, 

b'n = -ib'uf - ^{bu? + «r + 1^2)'. (Al) 
1 



ill — 



2b1iz 



(A2) 



+ 16<2<i: (A3) 
+2/f2&^2 + 2&?2/r2 - bl2fi2, (A4) 

f?2 = (b'^r + ^ibi2r + n<ir + «2r 

+ IK2)\ 

fi2 = 25?26?2 ~ i(6?2)2 - (/- )2 + 8u1,ul, 

-2K2)', 



(A5) 

(A6) 
(A7) 



'12 "12 1 (A8) 
<2 = -26nii?2 + 26?iu?2 - &ii"i2 + 46?2"n 

-f 2/- 7.^2 + 2/f2^?2 - fl2<2 ■ ( A9) 

Here g = 2TTvdg/dl with b = e'^' the dimensionless re- 
caling length of the RG transformation. The last three 
flow equations describe the renormalization of momen- 
tum non-conserving Umklapp processes. 

APPENDIX B: GAUGE REDUNDANCY 



operators were bosonized. Local gauge transformations, 
(^pia — *■ (f'pia + '^''^ Npia{x , t) for integer Npia leave the 
electron operators invariant, and so do not change the 
physical state. Thus, any shift in 9a and ipa which corre- 
sponds to an integer shift in (j)pia/2'K is redundant, and 
leaves the physical state unchanged. 

To establish whether or not two different semiclassical 
solutions, 9a and 9a, are actually physically equivalent 
we proceed as follows. For the given (spatially constant) 
shift S9a = {9a — 9'^)/2'K, we ask whether it is possible 
to choose appropriate (spatially constant) shifts Stpa so 
that the chiral flelds <j)pia/2'K are changed by integers. 
The choice for 5ipa is uncontrained, since the full inter- 
acting Hamiltonian is invariant under arbitrary spatially 
constant shifts in the four Lpa fields. If it is possible, then 
the two semiclassical solutions are physically equivalent. 
For physically inequivalent states, it will not be possible 
to choose Sip a to give the required integer shifts. 

To implement the above procedure we need an expres- 
sion relating the bare chiral fields (j)pia to 9a and Lpa ■ This 
can be obtained from 



1 



{ipp+ + atp. 



qaipa 



qa9^_ - q9p_), 



(Bl) 



where q = (—1)* = 1,-1 for bonding and anti-bonding 
bands respectively. For the D-Mott phase the relation 
between 9aJ£a and the fields 6'^±, iy9^± is given explic- 
itly in Eq. |3.3| . In the S-Mott phase, the eqivalence of 
(?4 = ipp- is modified to be 9f = (pp- + tt, but this tt 
difference does not effect the shifts 59a between different 
semiclassical states. Thus the ground state degeneracies 
in the D-Mott and S-Mott phases must necessarily be 
the same. In subsection B.l below we show that both of 
these phases have unique ground states. It is necessary to 
consider the SP and CDW phases separately (in subsec- 
tion B.2 below), since there is a non-trivial modification 
in the relation between 9a, fa and the fields 9p±, (p^±. 
As we shall show, in these latter two phases the ground 
state is two-fold degenerate - corresponding physically to 
the spontaneous breaking of a discrete parity symmetry 
(see Section VI). 



The Bosonized sine-Gordon form for the SO (8) Gross- 
Neveu model appears to have a highly degenerate ground 
state. In terms of the four Boson fields 9a, the 
semiclassical ground states correspond to spatially uni- 
form values chosen to minimize the potential V{9) — 
— 5 ^^_^{, cos(6'a) cos(0b). Solutions include 9a = 2nna 
as well a.s 9 = 2'Kna + tt for arbitrary integers Ua- But as 
we shall see, in most situations these multiple solutions 
actually correspond to the same physical state. To see 
which solutions are physically equivalent, it is necessary 
to relate the 9a and their dual fields (pa to the original 
Boson fields, 4>pia, introduced when the electron fermion 



1. D-Mott and S-Mott phases 

In the D-Mott and S-Mott phases, shifts in the fields 
9a and ipa induce shifts in the chiral fields, 5<j)pia, of the 
general form 

Hp^o. = \p{Ap)abS9, + ^{Ap)ab5ipb, (B2) 

where a=l|, 11,21,24 labels the band and spin in- 
dices and b — 1 , . . . , 4 labels the four flavors of the sine- 
Gordon bosons. Here and below, all shifts will be mea- 
sured in units of 2tt, so that for example 6(j) = (0— (/)')/27r. 
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The matrice s Ap can be explicitly constructed by using 
Eqs. O, 




(B3) 



It will be convenient to separate out the two contribu- 
tions coining from the shifts in 9a and (jfa, respectively, 
by defining 



6^, 



{Apr'^Sip,. (B4) 



Comparing two semiclassical solutions, 9a and 9a, de- 
termines the shifts 6Qpa- These two solutions are physi- 
cally equivalent, provided shifts S(pa can be chosen so that 
the following eight constraint equations are satisfied, 



epa+S^Pa=mpa, 



(B5) 



with integer Npa- In this case, all eight shifts 5(j)pia are 
integer, and the electron fields are left unchanged. 

Since the four shifts, 5tpa, determine both right and left 
vectors, 5$Ha,^^*i,a, these two vectors are not indepen- 
dent, and similarly for the 9 shifts. Indeed one can see 
that. 



1111 









or. 






iM)ab = 









(B6) 
(B7) 

(B8) 



We can now use the eight constraint equations to elim- 
inate (5$ and arrive at four equations for SO. To this end, 
upon multiplying by M on the left sector of Eq. B5, one 
obtaines (5$fja — SQ^a = 4(iW)a{,7Via. Upon combining 
with the right sector of Eq. MB one obtains. 



Se^a = 2Naa - {M)abN^b. 



(B9) 



Two semiclassical solutions (which determine (56 ^a) are 
then physically equivalent provided these four equations 
have solutions for integer Npa- 

All of the semiclassical solutions take the form 9a ~ 
2naTr or 9a = {2na + l)7r with arbitrary integers ria- It 
is straightforward to show that for any two of these so- 
lutions the difference 69a corresponds to S&aa which are 
either even integer for all a — I, ...,4 or all odd integers. 
When they are even integers, Eq. p39| can be solved for 
integer N^a by taking N^a = 0. For odd integer SQna a 
solution with integer Np^a is also possible by taking, for 
example, N^^a = Sai- 

We have thereby established the physical equivalence 
between all of the semiclassical solutions. This implies 
that the D-Mott and S-Mott ground states are unique. 



2. SP and CDW phases 

In the SP and CDW phases, the relation between 9a, fa 
and 9fj,±,iff^± are changed, so the above conclusions are 
modified. In particular, one has 



9r 



•fa-, 



(BIO) 



in the SP phase and 0f = 0f ^ + ^, ^f-uvv ^ qSP 
in the CDW phase. Because the Boson fields are defined 
differently, the matrix Ap which relates the two sets of 
fields in Eq. B2 is modified. The appropriate matrix in 



this case, denoted Ap, is given by, 



{Ap)ab = 



1 


1 


p 


p 


1 


-1 


-p 


p 


1 


1 


-p 


-p 


1 


-1 


p 


-p 



(Bll) 



Notice that Ap = Ap, although the left matrices dif- 
fer in the third column. Similarly, the matrix M is also 
modified, becoming 



{M)ab 



/O 2 

2 

2 

Vo 2 0, 



(B12) 



Physical equivalence between two semiclassical solu- 
tions for the SP or CDW phases is, once again, estab- 
lished by finding a solution of Eq. B9 with integer Npia , 
except with M replacing M . As before, the difference 
between any two of the semiclassical solutions leads to 
either even integer or odd integer 5Qpa. For even inte- 
ger 5Qpa a solution is again possible by taking N^^a — 
and choosing appropriate integers for Npa- However, a 
solution in the integers is not possible for two semiclas- 
sical solutions_ differing by an odd integer shift vector 
5Qpa (since MabNpb is always odd). Two such semi- 
classical solutions would thus correspond to physically 
distinct phases. 

The fact that the ground state is actually two-fold de- 
generate can be established as follows. Consider two spe- 
cific semiclassical solutions, 9\ = Q and 6'^ = ^irSai- One 
can readily show that the shift vector, 6Q\^, connecting 
these two states is an odd integer vector, so that these 
two states are physically distinct. Next consider an ar- 
bitrary third semiclassical solution, 9a- If the relative 
shift vector between the first and third solutions, SQ\'^ is 
even then the physical states are equivalent. If, on the 
other hand, SQa^ is an odd integer, then (58^'^ is neces- 
sarily even, and the second and third solutions describe 
the same physical state. It is thus clear that there are 
only two physically distinct ground states in the SP and 
CDW phases. As discussed in Section VII this two-fold 
degeneracy corresponds to a spontaneous breaking of a 
discrete parity symmetry. 
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APPENDIX C: GAMMA MATRICES AND 
SPINOR REPRESENTATIONS 

In this appendix, we discuss some technical details of 
gamma matrices and spinor representations of S0{5). 
In general, there are two types of representations for 
SO{N). The first are tensors, which transform like 
products of vectors. Irreducible representations are then 
found by taking symmetric and anti-symmetric combina- 
tions (Young tableaux). However, to describe how (com- 
plex) fermions transform under rotations, the second rep- 
resentation, the spinor one, is necessary. It has already 
been used in constructing invariants in Sec. VI, but here 
we review the mathematics in somewhat more technical 
detail, in order to allow the reader to perform concrete 
calculations if he or she so desires. To explain the spinor 
representation, let us introduce a set of N generalized 
Dirac matrices which obey the Clifford algebra, 



(Cl) 



where A, B = 1,2, . . . , N . We then construct the N{N - 
l)/2 generators defined as commutators between all pairs 
of these Dirac matrices. 



(C2) 



It is easy to show that these generators satisfy the SO{N) 
commutation relations 

(C3) 

For N = 5, we choose a specific set of matrices to repre- 
sent the 5*0(5) group. The minimum dimension of a set 
of five matrices which satisfy the Clifford algebra is 4 x 4. 
Our particular choice is 







1(7, 



a- 

a* 



(C4) 

A useful property of the spinor representation is its 
"reality" . This means that the conjugate representation 
_^pafc^* ^Yso obeys the algebra, and is equivalent under 
a unitary transformation to the original representation. 
This follows because it is always possible to find a matrix 
R which satisfies the properties 



-1, R-^ =R'^ = R^ ^ -R, 

-(r-^^)*, R-^v-R^ {V-y. (C5) 



For N with our particular choice of Dirac matrices 
in Eq. C4, the matrix R is simply 



R 



1 

-1 



(C6) 



where 1 is the two by two identity matrix. The matrix 
R is useful in constructing irreducible representations of 
50(5). 



As we have seen in Sec. VI, these abstract matrices can 



be elevated to physical operators by sandwiching them 
between two spinors. The useful details are already given 
in the text of Sec. VI. Here we provide some reasoning 
and motivation for the choice of spinor taken there. For 



convenience, we copy the spinor definition from Eq. 6.7 



*(fc) 



aa{k) 



N) 



(C7) 



where N = (tt, tt). Here (/)fc is a complex function with 
absolute value one, chosen by Rabello et. alo to have 
D-wave symmetry in two-dimensions. As discussed in 
Section VI, this factor plays no role in the case of the 
two-leg ladder, and can be set to unity. At first blush, 
the particular choice of spinor appears rather arbitrary. 
It is not, for several reasons. At half filling, the system is 
particle-hole symmetric. For every hole excitation at mo- 
mentum k created by a{k), there is a particle excitation 
counterpart at momentum k — N created by a^{k — N). 
Parity symmetry implies there is also a particle excita- 
tion at the opposite momentum —k + N. Because these 
excitations occur symmetrically, they are chosen as up- 
per and lower components in the spinor. The use of the 
parity symmetry is not essential. However, it is rather 
convenient for later analysis in weak coupling because, by 
such a construction, all components have the same chiral- 
ity, i.e. act on the same side of the Fermi surface. Since 
the four-components spinor "^{k) contains excitations at 
both k and —k + N, the momentum k is only allowed 
to run over half of the Brillouin zone. The halved region 
in momentum space is also known as the folded Brillouin 
zone (shown in Fig. 9). Finally, one would like the spinor 
to obey canonical anti-commutation relations so that it 
annihilates or creates fermionic excitations. This is the 
origin of the constraint on 0^ : direct calculation verifies 
that, provided \ (l>k\'^ — 1, the canonical anti-commutation 
relation is satisfied. 



{^,{k),^l{k')} = {27rr6ab5{k-k'). 



(C8) 



Further straightforward algebra demonstrates that when 
the spinor sati sfies cano nical anti-commutators, the cur- 
rents in Eqs. 6.14 - 6.18 satisfy appropriate Kac-Moody 
algebras. This exercise, which we do not reproduce here, 
verifies that these currents are indeed 50(5) scalars, vec- 
tors, and tensors, as indicated in Sec. [v| . 

We conclude this appendix by obtaining expressions 
which relate the 28 50(5) currents defined in Sec. VI, to 
the 28 50(8) currents, G^^, introduced in Sec. IV. These 
relations can be determined by bosonizing the 50(5) cur- 
rents, rewriting in terms of the GN bosons 0a and ipa 
and the Klein factors, refermionizing, and finally chang- 
ing from Dirac to Majorana fermions. For example, 
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(C9) 



The general relations can be conveniently presented in 
the following form: 



^ p 

Sp 



(CIO) 



AB 



The 5x5 anti-symmetric tensor matrix T^^ = for 
B and it is zero for A = B. The 3x5 vector matrix 

is 

v?- = i(-iSr-| . (cii) 

Finally, the 3x3 anti-symmetric scalar matrix is 

Sp = ^ \ Rclp ) . (C12) 
Plmlp PJp 



APPENDIX D: S'0(5) CURRENTS IN SU(2)xU(l) 
AND SO(8) NOTATION 

In Section VI the most general set of 5*0(5) invariant 
interactions for the weak coupling two-leg ladder were 
expressed as pro ducts of right and left moving 5*0(5) 
currents, see Eq. 3.20. Here we re-express these five in- 
teractions in terms of charge and spin currents with lower 
C/(l) X SU{2) symmetry, which were introduced in Sec- 
tion II. The products of 5*0(5) scalar, vector, and tensor 
currents are re-expressed as 



1^ f?. L 



{Jr.11 ~ Jr22 ~ 2)(Jj^ii — Ji22 — 2) 
4( J nil — J r22){J Lll — J L22) 
+ 2(^Il2-^i21 +-^121-^^12)1 



1 



Jr" Jl" — -^{JrII + '^r22 — 2)(Jill + Ji,22 — 2) 
+ 2 (J Rll + J r22){J Lll + J L22) 
-Mlil2h21+ll2llRl2)- 



(Dl) 
(D2) 

(D3) 



Notice that these three interactions conserve the num- 
ber of particles in each band. The remaining two 50(5) 
invariant interactions, involving anomalous scalars and 
vectors, scatter particles from one band to the other. In 
terms of the f/(l) x SU{2) charge and spin currents, they 
are 



I^I^ + Till = 4( Jp21 Jl21 + J«12 J^12), (D4) 

l^lt +If If = 16(Jr21 J.21 + JR12J.12) 

~2(/lii/j^22 + Jr22Jl11 
+4ll422+422^Rll)- (D5) 



For a given set of five 5*0(5) invariant interaction pa- 
rameters, these operator identities enable us to obtain 
the corresponding values of the nine forward, backward 
and Umklapp scattering amplitudes; 



b'li- 

K2 = 

fl2 = 
„P - 



4/1., 



5?2 = -16K„ 
fi2 ^ 45„ - 2gt, 



2K, 



,p _ 



9v, 



'■12 



(D6) 
(D7) 
(D8) 
2gt. (D9) 



From these, and the nine RG flow equations in Appendix 
A, one can obtain a closed set of five RG flow equations 
for the five 5*0(5) invariant coupling constants, given ex- 
plicitly in Appendix E. 

It is also instructive to re-express the five 50(5) invari- 
ant interactions in terms of the 50(8) currents - specifi- 
cally the 28 50(8) generators G"*^ = 177^77^, comprising 
the vector (fundamental) representation of 50(8). For 
the first three 50(5) interactions one finds. 



L 5 



5 

J'AB rrAB \ ^ r^AB /^AB 

R Jl — 2^ ^R 



(DIO) 
(Dll) 

(D12) 



A,B=1 



As expected, these expressions show that C"* , G"*** and 
QAB j-j-Qj, A, B = 1, ..,5) transform under 5*0(5) rotations 
as scalar, vector and (rank two) tensor, respectively. The 
remaining two anomalous 5*0(5) invariant interactions 
can similarly be re-expressed as, 



(D13) 



q-Aq-A , q-A'iq-A'i 



A=l 



{Gi'Gf + GfGf). (D14) 



It is clear that G"" '=** and G-*" 
scalars and vectors, respectively. 



transform as 50(5) 



APPENDIX E: 5*0(5) RG EQUATIONS 

For the weakly interacting two-leg ladder at half- 
filling, requiring 50(5) symmetry reduces the number 
of marginal four-fermion interactions from nine down to 
five. Due to symmetry, one expects the RG flow equa- 
tions to close in the manifold of 5*0(5) invariant models. 
This closure can be demonstrated explicitly by combin- 
ing the expressions obtained in Appendix D which spec- 
ify the five dimensional 50(5) invariant manifold, with 
the general RG flow equations in Appendix A. When re- 
expressed in terms of the 5*0(5) couphngs, the nine RG 
flow equations are seen to be redundant - only 5 are inde- 
pendent. Thus confined to the 5*0(5) invariant manifold, 
the set of independent RG flow equations can be written. 
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gv = S>gvgt - 52hshy, 

gt = 8gl + 6gf + 64hl, 
hs = ~^gshs ~ 2Qgyhy, 
hv = ~^gvhs - 'igshy + Sgthy. 
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